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At high energies the particles move very fast so their trajectories can be approxi- 
mated by straight lines collinear to their velocities. The proper degrees of freedom 
for the fast gluons moving along the straight lines are the Wilson-line operators 
- infinite gauge factors ordered along the straight line. I review the study of the 
high-energy scattering in terms of Wilson-line degrees of freedom. 
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1 Introduction 



Traditionally, high-energy scattering in perturbative QCD (pQCD) is stud- 
ied by direct summation of Feynman diagrams. In the leading logarithmic 
approximation (LLA) 

a s < 1, a s In — ~ 1 , (1) 

m 

the amplitudes at high energy are determined by the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) pomeronEl (for a review, see Ref. 2), 

/ a \ 12^ In 2 

" • (2) 

Here m is the characteristic mass or virtuality of scattered particles (for ex- 
ample, for the small- a; deep inelastic scattering m 2 = Q 2 ). In order for per- 
turbative QCD (pQCD) to be applicable, m must be sufficiently large so that 
a s (m) <C 1. 

The power behavior of BFKL cross section (0) violates the Froissart bound 
and, therefore, the BFKL pomeron describes only the pre-asymptotic behavior 
at intermediate energies when the cross sections are small in comparison to the 
geometric cross section 2irR 2 . In order to find the true high-energy asymp- 
totics by analysis of Feynman diagrams we should sum up not only the leading 
logarithms (a s Ins)™ but also the sub-leading ones a s (a s Ins)™, then the sub- 
sub-leading terms a 2 (a s lns)™, etc. This is almost equivalent to finding an 
exact answer to arbitrary QCD amplitude in all orders in perturbation theory. 
A more realistic approach is to unitarize the BFKL pomeron, i.e. to sum up the 
subset of sub-leading logarithms which restores the unitarity in s channel. Still, 
it is a difficult problem which has been in a need of a solution for more than 20 
years. One of the most popular ideas on solving this problem is reducing QCD 
at high energies to some sort of low-dimensional effective theory which will be 
simpler than original QCD, maybe even to the extent of exact solvability. The 
first step on this road is to identify proper degrees of freedom for this effective 
theory. One of the possible-choices is to formulate high-energy scattering in 
terms of "reggeized gluons.'o An alternative and related approach t£l is based 
on so-called Wilson lines - infinite gauge links corresponding to fast gluons 
moving along the straight-line classical trajectories. 

An important aspect of the Wilson-line approach to high-energy scatter- 
ing is the fact that it serves as a bridge between pQCD calculations and the 
semiclassical approach to high-energy scattering based on the solution of the 
classical equations for the fast-moving sourceso The semiclassical QCD (sQCD) 
is applicable when the coupling constant is small but the characteristic fields 



3 



produced by colliding particles are large, ~ ~. As advocated in Ref. 4, sQCD 
may be relevant for the heavy-ion collisions because the coupling constant can 
be relatively small due to high density of partons in the center of the collision. 
The relevant "saturation^ scale" was estimated to be ~ 1 GeV at RHIC and 
-2-3 GeV at LHC&BQ 

Let us demonstrate thatithe relevant degrees of freedom for the high-energy 
scattering are Wilson lines H As a result of the high-energy collision, we have 
a shower of produced particles in the range of rapidity between those of the 
colliding particles. Consider two clusters of particles with different rapidities: 
"A" particles with rapidities close to t\a and "B" particles with rapidities 
~ t\b- From the viewpoint of the "B" particles the "A" gluon moves very 
fast, so its trajectory can be approximated by a straight line collinear to the 
gluon momentum, see Fig. [|. The propagator of such gluon reduces to the free 

U(x , n A ) 

> 

' ' Wilsonline"-infinitegaugeliiik 



fast Ol~Ti A ) 




slow (rj-rj ) 



Figure 1: Propagator of a fast "A" gluon in the slow "B" background. 

propagator multiplied by the infinite gauge factor (made from "B" gluons) 
ordered along the straight line parallel to ua, the direction corresponding to 
the rapidity tja- 

U(x, ua) = [oon^i + x, — ooua + x]. (3) 
Hereafter we use the notation 

[x,y]=Pexpig du(x - y)" A^ux + (1 - v)y) (4) 
Jo 

for the straight-line gauge link connecting the points x and y. Therefore, the 
B particles can interact with A fields only via the Wilson lines (||) . Similarly, 
if we sit in the rest frame of the "A" gluons the "B" particles are moving fast 
along the direction collinear to the vector tib corresponding to rapidity t]b , see 
Fig. ||. The propagator of these gluons reduces to the Wilson line (made from 
"A" gluons) collinear to tib 

U(x, ub) = [ooub + x, — ooub + x]. (5) 

Again, the relevant degree of freedom is the non-local Wilson line @ rather 
than the local field A(x). We see that the particles with different rapidities 
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slow (T|~T| ) 

U(x a , n B ) 
fast (ri~ri B ) 

Figure 2: Gluon of "B" type viewed from the rest frame of "A" gluons. 

perceive each other as Wilson lines. The formal proof of this statement in 
terms of Feynman diagrams is given in the Appendix (see also Ref. 9). 

In this review I give a pedagogical introduction to the Wilson-line-based 
approach to high energy scattering. After a short overview of the traditional 
approach, I shall present the operator expansion for high-energy scattering 
which provides the operator language for the BFKL equation in the same way 
as the usual light-cone expansion gives the operator description of the DGLAP 
equation. Unlike the latter, there is a symmetry between the coefficient func- 
tions and matrix elements in the high-energy operator expansion which can 
be summarized by the factorization formula for high-energy scattering. This 
factorization formula gives us the rigorous definition of the effective action for 
a given interval of rapidity. In the last section we discuss the semiclassical 
approach to effective action related to the problem of scattering of two shock 
waves in QCD. 

2 The hard pomeron in pQCD 

Since there are many excellent reviews of the traditional, Feynman diagrams- 
based, approach to high-energy scattering (see e.g. Refs. 2, 10), I will present 
here the short introduction to the subject so as to set up the stage for the 
subsequent analysis of the high-energy scattering in terms of Wilson-line op- 
erators. 

2.1 High-energy 7*7* scattering 

For simplicity, we consider the classical example of high-energy scattering of 
virtual photons with virtualities ~ — m? 

A(s,t)=-i J d 4 xd 4 ^ 4 ze^^^^^^(0|r{^(x)j B ( 2/ )^(z)^(0)}|0). (6) 

Here ja{x) is electromagnetic current j^{x) multiplied by the polarization 
vector e^(p). In the Regge limit (s ^> m 2 , t) it is convenient to use the Sudakov 
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decomposition: 

P» = Optf + M+ft, (7) 
where p± and p^ are the light-like vectors close to pa and pb-, respectively: 

2 2 
T) T) 

PA=Pl + —P2, PB=P2 + —Pl, r = p B -p' B = a rPl + PrP2 + T±. (8) 

s s 

The momentum transfer r = p' A — pa = & r pi + PrP2 + fx has components 
a r r~j (3 r ~ — so t ~ — r . The typical diagram for the high-energy 7*7* 
amplitude is shown in Fig. 3 (recall that the diagrams with gluon exchanges 
dominate at high energies). 




Figure 3: A typical Feynman diagram for the high-energy 7*7* scattering. 

We will calculate the imaginary part of the amplitude A(s,i) 

W = -ImA. (9) 

The real part of A(s,t) can be restored using the dispersion relations. (It 
turns out that in the leading logarithmic approximation (LLA) the amplitude 
at high energy is purely imaginary, see e.g. the review in Ref. 2). 

Let us start with the lowest-order diagrams shown in Fig. 4. The integral 
over gluon momentum k = ctkPi + PkP2 + k± has the form 

W ° = ^ 4 / ^^^^M^k,H--k)lm^{-k,k-r) (10) 

where ($^)|^(fc, r — k) and (3?b)^( — k, k — r) are the upper and the lower 
blocks of the diagram in Fig. || (stripped of the strong coupling constant g). 
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J 



Figure 4: Lowest-order diagrams for the high-energy scattering of virtual photons. 



Here a, b and £, 77 are the color and Lorentz indices, respectively. In the Regge 
kinematics (= s ^> everything else) at ~ ^- and /?£ ~ x so fc 2 ~ — k\. 
Moreover, alpha's in the upper block are ~ 1 so one can drop in the upper 



block. Similarly, beta's in the lower block are 
the lower block. We get ($ afc = ^y$ cc ) 



1 hence one can neglect j3k in 



W u = 



2g 4 



(N? - 1)tt 
d 4 k 1 



(11) 



Im $^(fc, r - k)\ a =Q Im $^ hb (-fc, jfe - r) 



where 7V C = 3 is the number of colors. At high energies, the metric tensor g 11 " 
in the numerator of the Feynman-gauge gluon propagator reduces to ~p%p%, 
so the integral ( pi] ) for the imaginary part factorizes into a product of two 
"impact factors" integrated with two-dimensional propagators 



W° 



where 



N J -I 



/ B (-fc ± ,-r ± ) = . 



5 »j 
P2P2 



S (N? - 1) 

S v 
PlPi 

;(7V 2 -1) 



d 2 kj_ 1 



1 



4tt 2 k\ (r-k)l 



E«2 



2vr 



7 A (fc ± ,r ± )/ B (-fc ± ,-r ± ), 
(12) 

, (13) 



Im*S(Jfe,r-Jfe) 



a a 



(— k, k — r) 



,(14) 



0k=O 



and (J] e 2 ) is the sum of squared charges of active flavors. The photon impact 
factor is given by the two one-loop diagrams shown in Fig. ||. 
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^ 



A 



^ 




5 
k 



r-k + 

Figure 5: Photon impact factor. 

The standard calculation of these diagrams yields0 
I A (kx,rx) = I A (kj_,rj_) — I A (0, rj_), 




(15) 



where 
I A (k x ,r ± ) 



1 f 1 da f 1 da' 



2 ,/q 2tt Jq 2i 



f a' r ^ - -i-i 
- [Pla'a* - \p\5l + ( P ' A ) 2 a']aa) (16) 

x 2aa)(l - 2a'a')Pl{e A ,e' A ) A _ + 4:aaa'a'[Pl (e A ,e* A ) ± 

- 2(P,e A ) ± (P,e" A )i_] + (e A ,e' A )±(p A - (p' A ) 2 )aa(l - 2aa) 
x (1 - 2a') + 4afi(l - 2a)a'(P, e A )j_ (r, e' A )± } 



for the transverse polarizations A, A' — 1, 2. Here P± = k± — r±a and (a, b)± 
denotes the (positive) scalar product of transverse components of vectors a 
and b. 

2.2 The BFKL kernel 

In the next order in perturbation theory there are two types of diagrams for 
the 7*7* amplitude: diagrams with 5-particle cut describing the emission of 
an extra gluon and diagrams with 4-particle cut as in Fig. 4 but with an extra 
gluon loop. 

Let us at first consider the diagrams with the 5-particle cut shown in Fig. 6. 
The contribution of the diagram shown in Fig. 6a has the form 



W, 



(5) 

(a) 



= ~9 

7T 



16tH 



d 4 k f d 4 k' Im$7 ab (/fc, r-k) lm^ mn (-k', k' - r) 



£r/mn / 



16tt 4 k 2 (r-k) 2 



(k') 2 (r~k') 2 



x f amc f bnc T^(k, k')2w5((k - k') 2 )6(a k )r V7ia {r -k,r- k') 
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k-k 




(a) 



r-k 




p-k+k 




p-k+k' 



(d) 



+ k 




where 



p'-k+k' 




p'-k+k 




(f) 



Figure 6: The efEective vertex in LLA. 




(h) 



lm^' mn (-k',k' -r) 
{k') 2 {r-k') 2 : 



T^\(k, k') = (k + k')xg^ v + (k' - 2k) lJ g\ t _ t + (k — 2k') ll g v \ 



r-k 
r-k' 



(17) 



(18) 



is the three-gluon vertex divided by g. (Strictly speaking, in order to obtain 
<I>^ and $ B we must add the diagrams with permutations of the quark lines, as 
in Fig. 4). As mentioned above, it is convenient to use Sudakov variables (^): 
k = otkPi + (3kP2 + k±, k' — a' k pi + P' k p2 + k' ± . We will see that the logarithmic 
contribution comes from the region 



l>a»c/~— , _~ / 3« / 3'«l, kl~(k' ± ) 2 



(19) 



In this region k 2 = ak(3kS — k\ ~ — k\. In the same way, (k') 2 = -~(k' ± ) 2 , 
(r — k) 2 — — (r — and (r — fc') 2 = — (r — k')\. As we mentioned above, at 
high energies we can replace g^ v in gluon propagators connecting the clusters 
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w, 



with different rapidities by 2^p-. With these approximations, the integral 
( |l7| ) reduces to 

(5 ) = 2 6 /2\4 r da k dp k d 2 k ± r da' k dP' k d 2 k' ± 1 1 
W ~ n 9 \s) J 4tt 2 4tt 2 J 4tt 2 4tt 2 %\ (f- k)\ 

x Im$y Qb (fc, r - k)f amc f bnc T„ °{k, k')2n6(a k (3 k s + (k - k'f x )9(a k ) 
x T.*Jr-k,r-k')^ X - — Im$^ m "(-fc',fc'-r). (20) 

Since a in the upper block is ~ 1, one can neglect a^-dependence in &a which 
leads to the replacement of J d(3 k ^A by the impact factor I (k±,r±), see 
Eq. (|l|). Likewise, / da k § B — > I B {k'±,r±) so we get 

r(B) _ 2g 6 N C (N 2 - 1) 2 n 2 



= v^^l«J (21) 

, j- I A (fcx,rj.)^ ^ =, -^I B (k' ± ,r ± ) 

^^r #st fc')2^K/3^s + (fc - #)i)0(a fc )r.. ;ff (r - k,r - k'). 

Let us now turn to the diagram shown in Fig. 6b. Since the gluon with mo- 
mentum k—k' now connects parts of the diagrams with different rapidities, we 
can replace g^" in this propagator by 2 I>2 ^ 1 . After that, the quark propagator 
with the momentum p + k' in the upper block reduces to 

f* - K l a ^i|^t^ 2 r W^-J-rr, (22) 

-(op + <)(/3p + /%)s + (p - fc')i - le -P h - ^ 

(recall that a p ~ 1, fi p ~ — ) . We see that in the transverse space this propa- 
gator shrinks to a point so the answer for the upper block is again I A multiplied 
by „, -, - . (The eikonal factor al \ . is the Fourier transform of the first term of 

the expansion of Wilson-line propagator (g) in powers of "external slow field" 
represented by gluon with momentum k'). The right part of the diagram in 
Fig. 6b is identical to that in Fig. 6a so we obtain 

Pk^ 



wffi = ^{^/^/^^(^■^ (23) 



1 27rJ(a fc /3^ + (fc - fc')l)^K)r.*.(r - fc, r - fc') 
' ' / B (fcl,r ± ). 
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The contribution of the diagram in Fig. 6c is calculated in a similar way. One 
can replace 

t c i>2 [K - a k + oQh + i>± - V± + IfiAt 9 ] , tc ^ i t a (24) 



(op -a k + a' k ){f3 p -(3 k + j3' k )s -(p-k + fc')l + »c $ _ ie 

and, therefore, 

. - ^ / ^fg^ (t± ,,. ±)/ - (25) 

-27r*(a fc) 9 / fc a + (fc-fc / )iMttfc)r.*.(r-fc,r-fc / ) 



& (r - fc) 2 



_i_ 

1 -B^ 



(fc')i (r - 



I"(k' ± ,r x ). 



Note that the sum of the results (^l|), j2^), and (|25| ) may be obtained from 
the contribution ( pl| ) of the diagram in Fig. 6a. by the replacement 

r./(fc,fc')^r./(fc,fc')-|Lp 2 CT . ( 26) 

Now consider now the the diagram in Fig. 6d. The two quark propagators 
carrying the momentum k' give 

(1 -a p + a' k )^ -j ± + lt' ± 

F2 



X H 



(1 - a p + a' k )(?f -(3 P + (3> k )s - (p- k')\ + ie 

(a P + a k )A + P'± + P± j 
(a p + a' k )(j3 p + (3' k )s- (p-k>)l+ie 
(1 - a p )i> 1 - fS ± + }6' ± A (a p + a' k )i >1 + tf ± + 



(1 - a p )f3 k s - a p/3' k s 



(27) 



Since we cannot keep both large terms (1 — a p )p\ and a p p\ in the numerators 
this expression is -S^- times smaller than the contribution (|2^) of the diagram 
in Fig. 6b so it vanishes in the LLA. 

The diagrams in Fig 6e,f are calculated in the same way as the diagrams 
in Fig 6b, c. Similarly, the result may be obtained from Eq. ( |2"T| ) by the re- 
placement 

r./fo*') (28) 
11 



In conclusion, the diagram in Fig. 6g vanishes in the LLA for the same reasons 
as the Fig. 6c diagram. 

Thus, the contribution of the diagrams in Fig. 6a-6e can be represented 
by one diagram shown in Fig. 6h: 



(5) sg 6 N c (N 2 c 

VV (a+...g) n 4 



-(E e ') 2 (»> 



J 4n 2 4ir 2 V x ' ± 'k 2 (r-k) 2 (k') 2 (r - k') 2 



(f-fc)i (fc')l(r- k') 2 ± 
^^L°(k, k')2TrS(a k p' k s + (k - k') 2 ± )9(a k )T.^(r -k,r- k'), 



where 



L"{k,k') = lTZAk,k')-2^p° 2 -2^pl 

s fi' k s a k s 

= {k + k')l-{a k + 2^-)pl-{(3' k +2 { ^)pl (30) 
P k s a k s 

is the Lipatov effective vertex for the gluon emission shown in Fig. 6h by 
a shaded circle. Note that unlike the usual three-gluon vertex, the effective 
vertex is gauge- invariant, 

(k- k') a L a (k,k') = 0. (31) 



We have demonstrated that if we take the diagram in Fig. 6a and attach 
the left end of the k — k' gluon line in all possible ways, the left three-gluon 
vertex in Fig. 6a is replaced by the effective vertex (^0|). Likewise, the sum of all 
possible attachments of the right end of this k — kl gluon line converts the right 
three-gluon vertex T.^ (r — k,r — k') into the effective vertex L a (r — k,r — fc') . 
Hence the sum of all the diagrams with 5-particle cut takes the form (see Fig. 7) 



W {5) 



4 U^J 2 



(32) 



' d 2 k ± d 2 k' ± I A (k ± ,r ± )I B (k' ± ,r ± ) 



Mi 



4tt 2 4tt 2 k2^_k)l(k') 2 x {T f -k')_ 
^^L°(k, k')27r6{a k 0' k s + (k- k') 2 ± )e(a k )L a (r - k,r - k'). 
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-X- 



k i i r-k 



Figure 7: Sum of the diagrams with gluon emission in LLA. Shaded circle denotes the 
effective vertex. 



Since a k (3' k s = — [k — k')j_ due to the 5-function, the product of two Lipatov's 
vertices gives 

\L'% k!)L a {r -k,r-k') = -ri + ^f-^f + ^r^ ' ^ 
1 [k — k')j_ (k — k')j_ 

which is proportional to the "emission" part of the BFKL kernel, see the 
Eq. (36) below. Now one can easily perform the remaining integrations over 
a k and j3' k in the LLA 

s [ da k d/3' k S{a k (3' k s + (k - k') 2 x )6{a k ) = [ da k — =ln-^, (34) 
J Jw2. a k m z 

and, therefore, the final result (for the diagrams with 5-particle cut) is 

WW = V^f^cln^ (35) 
7r 4 2ir ui 

00I A (*x,rx)^^Xx(*x,*l,r)/-(*l,rx) 
where 

k\ (k-f)l 



K(i)(k±,k' x ,r) 



(fc')i(r - k' ± ) 2 (k')l(k-k' ± ) 2 (fc' - r)l(k - k' ± ) 2 

(36) 

is the first part of the BFKL kernel coming from the diagrams with gluon 
emission. 
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(a) 



(b) 



(C) 



(d) 



+ k' 



n'-k+k' 



r-k + 



(e) 



p'-k+k 




r-k + 



p'-k+k' 



r-k 



(f) (g) 
Figure 8: Virtual corrections. 



k' k-k' i r-k 



(h) 



Apart from the diagrams with 5-particle cut shown in Fig. 6, there are 
also diagrams with four-particle cut ( "virtual corrections" ) of the type shown 
in Fig. 8. Let us consider the diagram shown in Fig. 8a. The integrals over 
aifc and [3k are similar to the same integrals in the first-order diagram in Fig. 4 

2 

and therefore ak ~ Pk ~ — • The logarithmic contribution comes from the 
region 1 ~ a p ^> a' k » ak- In this region we can replace the quark propagator 
with momentum p — k' by the eikonal propagator (see Appendix 7.1), 



i 



-P' k + ie- 



(37) 



In addition, one can neglect j3' k in comparison to (3' p ~ 1 in the lower block . 
The loop integral over k' turns into 



da' k d(3' k d 2 k' r 



47T 2 47T 2 



t a - 



a' k /3' k s-kla' k (0' k -p k )-(k-k')l 



(38) 
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The integral over j3' k is determined by the residue at (3' k — so we obtain 



2l 2na' 



d 2 k' ± 

1^2" 



t a fat 1 



1 



1 



k\ (k - k/) 2 ± 



t a fat b t a fat 



4tt 2 



In A 



d 2 k' ± 1 



t a fat b 
i 



™ 2 J 4tt 2 (fc-fc') 2 ^ 



(39) 



Let us add now the contribution of the diagram in Fig. 8b. Like the Fig. 8a 
case, we get the loop integral over k' in the form 



da' k df3' k d 2 k' 

47T 2 47T 2 



fa 



-P'k + ^ 



1 



1 



*' k P' k s-kl a> k ((3> k -(3 k )-(k-k>) 



a' k P' p3 - (f + k')\ 



da 'k 

d 2-rroi 



d 2 k' ± 
4tt 2 



t a fat 1 



k\ (k-k') 2 ± 



t a fat b t a faf 



Air 2 



In A 



d 2 k' ± 



t a fat b 
1 



4tt 2 



(k')l(k-k') 2 ± 



(40) 



The diagrams shown in Fig. 8c-g do not give the logarithmic contribution for 
the same reason as the diagram in Fig. 6d. 

We see that the sum of diagrams in Fig. 8a-g reduces to the first-order 
diagram in Fig. [|a with the left gluon propagator — p- replaced by the factor 



1 



-hi 



d 2 k' ± 



1 



4ir 

shown schematically in Fig. 8h. We get 

N 2 - 1 -n 2 - 2 



(fc')i(fc-fc') 2 



(4) 

(a+...g) 



= 9 

IT 



V 4tt m 2 



(41) 



(42) 



d 2 k ± I A (k ± ,r ± )I B {k ± ,r ± ) J f d 2 k' ± 



An 2 



kl(r-k)l 



(fc')i(fc- w) 



i\2 



The diagrams with the gluon loop to the right of the cut lead to similar re- 



placement of the right gluon propagator 



by 



-{k-f)l 



9 Z N C s 

~A 111 ~~0 

47T m- 2 



d 2 k' ± 



4tt 2 



(k' ± -r ± ) 2 (k-k f ) 



'\2 



(43) 
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Thus wc obtain the result 



^ (4) = -V^(£eS) a £jVoln4 



7T 4 V 47T m 2 

d 2 k ± I A (k ± ,r ± )I B (k ± ,r ± ) 
4tt 2 fc2( F _^ 

d 2 k' ± j P ^ jk-f) 1 

4tt 2 



(fc')i(fc-A:')i (fc'-r)2_(fc-fc') 2 



(44) 



for the contribution of the diagrams with 4-particle cut. 

Adding the sum of the diagrams with real gluon emission W 5 we obtain 
the final result for the 7*7* scattering amplitude in the first order in LLA. It 
can be represented in the form 



txi-1 3 4 N c- 1 /v- 2 \ 2 9 » r , a d 2 kd 2 k' 

W 1 = -a — > e 2 n ) —N c \n— — - — - (45) 

7r 4 V 2tt c m 2 / 4tt 2 4tt 2 V ; 



where 



N 2 ~ 1 ,\ 2 a 2 „ , s fd 2 kd 2 k' 
e 

x ^(^^x) 1 iC(fc ± ,fc^,r)J B (fc^,rx) 
/c_!_(r — fc)j_ 



K{k ± ,k' ± ,r) = K (1) (k ± ,tf ± ,r) - U^(k-k') (46) 

d 2 fc"_L (fc-r) 2 1 



is the BFKL kernel^ The explicit form of K is 

k\ (f-k)l 



K(k±,k\,r) = 7> 1 t = = =■. — 



^(f-ifcoi ^(fc-fcoi (?=*- fc')i( fc - fc ') 



_w- fc0 /^ , g2 , + , ( fc -y, . (47) 

2 7 4tt 2 \ - (fc" -f)i(A;-fe")iJ 

Note that both VF^ 5 ) and 1F< 4 ) are IR divergent but their sum W given by 
Eq. (45) is IR finite. This is the usual Bloch-Nordsieck cancellation between 
th emission of real gluon in diagrams in Fig. 6 and virtual gluon in Fig. 8. 
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2.3 Bare pomeron in the LLA 

The 7*7* amplitude in the first two orders in perturbation theory may be 
represented in the operator form as 

^(o +1) = sC r *k jA{ ±) i / t + l_ NcXn i_ k \jB {k ^ 

J 4tt 2 k\{r-k)\^ 8tH m 2 J 

(48) 

where C = a s (N^ — l)(J2 e q) 2 an d the operator K r is defined by its kernel 
K(k,k',r), 

(K r f)(k ± ) = J l^K(k ± ,k' ± ,r)f(k' ± ). (49) 

We can demonstrate (and we will do this using the evolution equations for 
the Wilson-line operators) that in the next orders in LLA the operator K 
exponentiates: 



It is convenient to represent the amplitude as an integral over the complex 
momenta: 

/•(5+ioo 



,0+100 w 

^(M) = ^-/ 5 _ ioo ^) ( 51 ) 

where u> = j — 1. The relation between the LLA and the power series for 
W{u,t) is 



W^(s,t) = S C^l(.g 2 ln-^)"/„(t) 

n—1 

00 2n 

V^ LA (c,t) = CJ2 /«(*) ( 52 ) 



where 



AW ^ /g'W V(*..rx) (53) 



are the coefficients of the LLA expansion. 
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The asymptotics of the amplitude at s — > oo is given by the rightmost 
singularity of the integrand in the right-hand side of Eq. ( |5l"| ) in the uj plane. 
The position of this singularity is given by the maximal eigenvalue of the 
operator K r determined by the eigenfunction equation 

^(K r f)(k ± )=u>f(k ± ). (54) 

This equation is solved at arbitrary momentum transfer r yet it turns out 
that the maximal eigenvalue of Eq. (50) does not actually depend on r. For 
simplicity, let us consider the case r = corresponding to total cross section of 
7*7* scattering. (In the next section we prove that the position of singularity 
does not depend on t = — r^). 

At r = 0, the full and orthogonal set of eigenfunctions of the BFKL oper- 
ator are simple powers 

f(k) = (fc 2 )^ + "V^, (55) 

with the eigenvalues 

a \n\ + 1 

u = 2Nc—x{v,n), x(v,n) = -Bs&( LL z—+w)-C. (56) 

7T Z 

The maximal eigenvalue is 2 J ?V c ^x(0,0) = 4 2jL iV c ln2, so the rightmost singu- 
larity (intercept of the "hard pomeron" ) is located at 

j = l+u , ^ = 4^iV c ln2, (57) 

so the asymptotics at high energies in the LLA is 

*-N c In 2 



(m 2 ) 



(58) 



It is easy to see that the singularity at to = ojq is the branch point -j==. 

As we mentioned in the introduction, the singularity at j > 1 violates the 
Froissart bound a < In 2 s. Recently, the next-to-leading correction (~ a s ) to 
the BFKL kernel was found^J but the result still violates the Froissart bound, 
so the unitarization of the BFKL pomeron is required. (Consequently, the 
BFKL pomeron ( |57j ) is sometimes called "the bare pomeron in pQCD" ) . 

In the case of 7*7* scattering, it is possible to find the explicit form of the 
cross section in the LLA. Expanding impact factors l(k,0) = I(k) in a set of 
eigenfunctions (|55|), we obtain 

a tot (PA,PB) = 5 4 i(iV c 2 -l)(£ e ') 2 (59) 
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Here we neglected the angle-dependent contributions coming from n =f since 
they decrease with energy. At s — > oo the cross section ( |59| ) is determined by 
the rightmost singularity in the v plane located at v = (in terms of j-plane 
it corresponds to Eq. (p7|)) and the result is 

where C(3) ~ 1.202. 



U.^ Diffusion in the transverse momentum and the BFKL equation with run- 
ning coupling constant 

At first, let us demonstrate that the rightmost singularity of the BFKL equa- 
tion is located at cu = ujq at t ^ as well (although its character changes 
from to y/uj — luq). We shall see that in higher orders in perturbation 

theory there is a "diffusion" in k± such that In ~ y/n (where n is the or- 
der of perturbation theory). To illustrate the diffusion, consider a rung of the 
BFKL ladder located in the middle of the rapidity region (see Fig. ^) . Each of 




Figure 9: Diffusion in fej_. 

the upper or lower blocks in this diagram are "non-integrated gluon distribu- 
tion" . The s — > oo asymptotics is governed by the rightmost singularity of the 
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function W(u),t) (see Eq. ([5l|)) which is determined by the asymptotics of the 
coefficients /„ at n —* oo. For even n, these coefficients can be represented as 

Mt) = J ^^^f-(k x ,r x )f-(k u r x ), (61) 

where 

f£(kx,r x ) = f^K r yi A (k x ,r x ), f*(k x ,r x ) = ^K r J I B (k x ,r x ). 

(62) 

Let us demonstrate that the characteristic momenta k\ in the integral in 

Eq. ( |6l| ) are ~ m 2 e v/ ™. At large transverse momenta k x the recursion for- 

A 

1 + 1 



mula fn+i{k X ,r x ) = ^±K r f^(k x ,r x ) can be reduced to 



w&H-i(0 = (63) 



(:' 



(?-«')/2 



1 - e«-f ' r v w V 1 ^ e? ~ ? ' Vl + 4e 2 («-« 



where £ = In ^ and </>«(£) = f V JF[J^ n (^-)" Next, we expand the function 
(j> n (C) m the integrand in Eq. ( p3| ) in Taylor series </>«(£') — <An(£) + {C ~ 
QftniQ + §(£' — £) 2 0"»(£) + As we shall see below, at large n and fcj_ one 
can neglect higher terms in Taylor expansion, and then the recursion integral 
equation ( |S3"|) can be approximated by the differential equation 

"iL^' €) = (W ° " W)0(n ' C) + c525 ^( n ' 0» (64) 

where c = -t-<7 2 C(3), C(3) — 1.202. This equation describes the diffusion of the 
"particle" where n serves as a time and £ as a coordinate. It is well known that 
at large time n the mean position £ of the "particle" is proportional to ^/ri, 
and therefore our approximation of Eq. ( |63| ) by the diffusion equation (|^) is 
justified. 

Thus, we must find the solution of the diffusion equation ( |64| ) with the 
"wall-type" boundary condition 

0(»,OI{= Ct =O, ^= ln ^ ( 65 ) 

which reflects the fact that our approximation is not valid at k\ < r^_. It is 
easy to check that the solution of the Eq. ( pi|) with the boundary condition 
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(pBI) behaves at large £ ~ sfn as 



^-^) e (?-l)n e -^(^-gt) 



where the coefficient of the proportionality may be determined by a more ac- 
curate analysis of the transition from the integral equation d63| ) to the diffusion 
equation (p4h. 

Substituting the estimate (pq) in the integral (pl|), we obtain 



jj 



1 



,3/2 



(67) 



which gives 



(2 \ /*oo i 

9 -\ f n = J dn-^e^- 1 ^ = VW^- (68) 

We see that the singularity is located at the same point u — luq as in the 
case of forward scattering, although its character is slightly different: ^ojq — oj 
instead of , 1 El 

At t = there is no "wall" boundary condition (pq ) which shows that the 
diffusion equation (|64|) leads to |£| ~ y/n. This means that the characteristic 
momenta k± are either very large, k\ ~ m 2 e v/ ™, or very small, k\ ~ m 2 e _v/ ™. 
The large contribution from the region of small k± region indicates the possi- 
bility of the breakdown of perturbative QCD for high-energy scattering. 

We can safely apply pQCD to high-energy scattering if the characteristic 
transverse momenta of the gluons k± in the ladder are large. For the 7*7* 
with p\ ~ p 2 A ~ m 2 ^S> Aq CD one can check by explicit calculation that the 
characteristic k± for the first few diagrams are ~ m. However, due to the 
diffusion in k± , the leading contribution to the loop integrals comes from 
the gluon momenta which are either very large, k 2 ^ ~ m 2 e v/ ™, or very small, 
k\ ~ m 2 e~^™. Due to the asymptotic freedom, the fact that the k± may 
be very large at n — » 00 only strengthens the applicability of pQCD. On the 
contrary, the fact that k± may be small questions the applicability of pQCD 
to the high-energy 7*7* scattering. 

To take into account the asymptotic freedom, one may consider the BFKL 
equation with the running coupling constant. Each of the upper or lower blocks 
in the diagram in Fig. |9| is a "non-integrated gluon distribution" 

F« B \k ± , r ± ; S )=E^ (a" ^ 4)" tf {B) (kx,r X ) (69) 
z — ' n\ V m z I 
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which satisfies the BFKL equation 

LjF A ( B \k ± ,r ± ;cj)= (70) 
I^ B \k^ ri _) + £^N C J d 2 k' ± K(k ± ,k' ± ,r ± )F A ^(k' ± ,r ± ;uj) 

where F(kj_,rj_;uj) is a Mellin transform of Eq. (]69|): 

F(k u r ± ;s) = ^ ( duj(AX ' F(k±,r±;u). 
2m J \m z / 

In order to account for the asymptotic freedom, we can replace g 2 in the right- 
hand side of the Eq. @ by g 2 (k 2 ± )^ 

u>F(k x ,r ± ;u) = I(fci,r x ) + l ^ N o J d 2 k' ± K(k ± , k' ± ,r x )F(k' ± ,r x ;u>). 

(71) 

This equation exceeds the LLA accuracy but it it can be demonstrated that 
in the case of large (or small) k\ the replacement g 2 — + g 2 (k 2 _) agrees with 
the renormalization group analysis E3). Another arguments in favor of taking 
into account these particular sub-leading logs follows from the analysis of the 
renormalon contributionsli-3 

At large k± one can replace the equation ( |7l| ) by the corresponding dif- 
fusion equation. It turns out that at large momentum transfer \t\ = fj_ the 
rightmost singularity of F(k±, rj_; uj) is located simply at t = 12 °*^*^ jV e ln2. 
At t = the diffusion goes in both directions leading to the contributions 
coming from k± ~ Aqcd ■ If one removes these contributions "by hand" (im- 
posing the "wall" condition at k\ — Aqcd), one obtains^ discrete set of Regge 
poles which condense from the right to the point u = OEj A more satisfactory 
solution of the problem of the diffusion to small k± would be to match the 
hard pomeron with the soft Landshoff-Donnachic pomeron (responsible for 
the high-energy hadron-hadron scattering) which presumably comes from the 
high-energy exchanges by soft gluons (see, however, Ref. 15 for an alternative 
"hard" soft pomeron). Another possibility is that the diffusion to small k± 
disappears if one takes into account the unitarization effectsll3 

The proper way to address the problem of running coupling constant in the 
BFKL equation is to use the NLO BFKL kernel in the renormalization-group 
analysisBJ The NLO correction to the anomalous dimension of the correspond- 
ing leading-twist gluon operator consists of two parts: the conformal part and 

'We have seen from the diffusion equation that (k' ± ) 2 ~ fc^ in the adjacent rungs of the 
ladder so g 2 (k\ ) = g 2 ((k') 2 ± ). 
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the running coupling part. The conformal part (see also Ref. 18) corrects the 

intercept of the BFKL pomeron (pfl), while the running coupling part, be- 

— 1 2 \ 

sides replacing 12^7V c ln2 by 12 ° 3 y ' N c \n2 in the leading order, leads to 
the non-Regge terms in the energy dependence of the cross section. The nu- 
merical value of the correction to the hard pomeron's intercept introduced by 
the conformal part of the NLO BFKL kernel is. large and negative. Its exact 
contribution is somewhat difficult to estimateuScd There are hopes, however, 
that collinear singularities causing this large NLO correction cancel each other 
at higher orders in a s c3 



2.5 Reggeized gluons and unitarization of the pomeron 

As I mentioned above, the bare pomeron violates the Froissart bound so we 
need to unitarize the BFKL pomeron. There are several approaches to the 
unitarization: effective reggeon field theorypj the generalized LLAcj eqiiijszalent 
to the quantum mechanics of reggeized gluons,^] and the dipole modelESEl We 
postpone the discussion of the dipole model until the next section and turn 
the attention to reggeon-based schemes of the unitarization. 

The reggeized gluon can be defined as a "hard pomeron" for the quark- 
quark scattering. We have seen that the gluon propagator describing the 

k x 

exchange between two quarks to the left of the cut in Fig. 4 is replaced in the 



next order by the factor (41) coming from two diagrams in Fig. 8a, b. Thus, 
in the first two orders in perturbation theory the propagator describing the 
exchange between two quarks with gluon (color octet) quantum numbers in 
the t channel has the form 

a s N c \n-^ I ? - ) . (72) 

m? J 4tt 2 (#)2(jfe-jfc')2j V ' 

It can be demonstrated (either by direct summation of the Feynman diagramsB 
or by evolution of the Wilson-line operators, see Sec. 3 below), that in the LLA 
the logarithmic factor in parenthesis exponentiates, therefore the exchange 
between two quarks is described by the "reggeized" gluon propagator 




(73) 



c In the reggeon quantum mechanics, the unitarity is preserved only in the direct s-channel, 
while in a reggeon field theory the unitarity holds true in all the sub-channels corresponding 
to different groups of particles in the final state. 
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where 



-a s N c 



d 2 k' ± 



(k')l(k-k'f 



(74) 



is the trajectory of the reggeized gluon in the plane of complex momenta in 
the leading order in a s ^\ Recently, this trajectory was computed in_the next- 
to-leading order in a s by direct summation of Feynman diagrams o and by 
calculatioiLof the two-loop anomalous dimensions of the relevant Wilson-line 
operatorsE3 

In terms of the reggeized gluons the BFKL ladder can be resummed as 
shown in Fig. 9 where the dash-dotted line denotes reggeized gluon ( f73f ) and 
the reggeon-reggeon-particle interaction is described by Lipatov's vertex (|30|). 
(The expansion of the reggeon trajectory in powers of g 2 reproduces the BFKL 



x=r) 



¥(p.,V 



Figure 10: BFKL ladder as a propagator of the two-reggeon state. Reggeized gluons are 
represented by dash-dot-dot lines. 

result ( |50| ) after combining the terms with like powers of g 2 ). This diagram can 
be interpreted as an evolution with respect to "time" = rapidity of the two- 
particle state described by the wave function ^(pi, P2) in quantum mechanics 

d This trajectory is IR divergent as it should be for the amplitude of the scattering of the 
colored objects. For the scattering of white objects (like virtual photons discussed in the 
previous section) this divergence will cancel with the IR divergence for real gluon emissions. 
To avoid the infinities in the intermediate results, one can use the dimensional regularization 
(with d = 2 + e transverse dimensions) or assume a small gluon mass fi. 
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with the Hamiltonian @ 

9 2 N C 



H 



12 



{ln|pi| 2 + ln|p 2 | 2 (75) 

+ J^-^(P*1P2^\P12\ 2 (P1P*2+C.C.)+4C} 

where pj = x±\ + ix± 2 i Pj — i-§f (index j = 1,2 numbers the particles), 
and pi2 is the coordinate operator (pi2 = pi — p 2 )- The first two "kinetic 
terms" correspond to the propagators of the reggeized gluons and the third 
term describes the interaction of reggeized gluons by exchange potential coming 
from product of two Lipatov's vertices given bjz_.Eq. ([36]). The Hamiltonian 
( [75| ) has a property of holomorphic separability^ 

H 12 = h 12 + h* 12 , (76) 



where 



hu = ttt-t |lnpij3 2 + (In p\2)P\ + — (In p 12 )$2 +2C\ , (77) 

167T 1 Pi P2 



and C=0.557 is Euler's constant. The generalized LLA is the summation of 
the diagrams shown in Fig. 11 (see the discussion in Ref. 29). The number 




*(P..tU 



X=T| 



Pl P 2 *(P.,V P n 



Figure 11: Generalized LLA as quantum mechanics of the reggeized gluons. 

of reggeized gluons in t channel is conserved, so the sum of the diagrams in 
Fig. 10a can be described by quantum mechanics of the reggeized gluons with 
pairwise interaction ([75|), 

H = Y,T?T£H ik , (78) 

i<k 



25 



where Hik is obtained from Eq. ( |75[ ) by the trivial replacement 1 — ► i, 2 — * fc. 

The unitarity follows from the representation of the sum of these di- 
agrams as a generalized eikonal E3 (see Fig. 12). In the multi-color limit 



exp <; 



Figure 12: Quantum mechanics of the reggeized gluons as a generalized eikonal. 



(N c — > oo, g 2 iV c -fixed), the non-planar diagrams vanish hence only the inter- 



action between the adjacent reggeons survives (the unitarity still hi 
The color structure is then unique and the Hamiltonian reduces to f 



H 



ds true). 



(79) 



where \ comes from the fact that the adjacent gluons are in the octet state. 
Using the property of the holomorphic separability ( |76| ) , it is possible to reduce 
the quantum mechanics of the reggeons described by the Hamiltonian ( |79| ) to 
the XXX Heisenberg model with spin s = OEil Unfortunately, the explicit 
solution for the number of the magnets k > 3 (= number of the reggeons) 
has not yet been found. For the k = 3 (the so-called Odderon state of three 
reggeized gluonsl±hfi variational estimates give the intercept at the value of J 
slightly below lE3o (recently, another Odderon- type solution with intercept 
at j — 1 was found in Ref. 34). 

In synopsis, we have found the subset of the non-LLA diagrams which 
restores unitarity in the s-channel and in the large N c limit this subset reduces 
to the one-dimensional quantum mechanical model (XXX magnet with s = 0) . 



3 Operator expansion for high-energy scattering 

The expansion of the amplitudes at high energy in Wilson-line operators is 
very useful in a situation like small-x DIS from the nucleon or nucleus. As 
the usual light-cone expansion provides the operator language for the DGLAP 
evolution, the high-energy OPE gives us the operator form of the BFKL equa- 
tion. In the case of deep inelastic scattering there are two different scales of 
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transverse momentum k±, and therefore it is natural to factorize the ampli- 
tude in the product of contributions of hard and soft parts coming from the 
regions of small and large transverse momenta, respectively. Technically we 
choose the factorization scale Q > fi > mpf, and the integrals over k\ > /i 2 
give the coefficient functions in front of light-cone operators while the contri- 
butions from k\ < fj 2 give matrix elements of these operators normalized at 
the normalization point \x. In the final result for the structure functions the 
dependence on \i in the coefficient functions and in the matrix elements cancels 
out yielding the Q 2 behavior of structure functions of DIS. 

In the case of the high-energy (Regge ) limit, all the transverse momenta 
are of the same order of magnitude, but colliding particles strongly differ in 
rapidity, thus it is natural to factorize in the rapidity space. Factorization in 
rapidity space means that a high-energy scattering amplitude can be repre- 
sented as a convolution of contributions due to "fast" and "slow" fields. To 
be precise, we choose a certain rapidity rjo to be a "rapidity divide" and we 
call fields with r\ > rjo fast and fields with rj < i]q slow where ryo lies in the 
region between spectator rapidity t\a and target rapidity tjb- (The interpre- 
tation of these fields as fast and slow is literally true only in the rest frame of 
the target but we will use this terminology for any frame). Similarly to the 
case of usual OPE, the integrals over fast fields give the coefficient functions 
in front of the relevant (Wilson-line) operators while the integrals over slow 
fields form matrix elements of the operators. For a 2=^2 particle scattering in 
Regge limit s 3> m 2 (where m is a common mass scale for all other momenta 
in the problem t ~ p\ ~ {p'a) 2 ~ Pb ~ (p'b) 2 ~ m2 ) this operator expansion 
has the formEJ 

Mpa,Pb p'aiPb) = ^ J d 2 x 1 ...d 2 x n C %1 - %n {x 1 , ...x n ) 

x (pB\TT{U ll {x 1 )...U tn (x n )}\p , B ). (80) 

(As usual, s = (pa +Pb) 2 and t = (pa ~p'a) 2 )- Here Xi (i = 1, 2) are the trans- 
verse coordinates (orthogonal to bothp^ andps) and Ui(x) = W(x)^-^-U(x) 
where the Wilson-line operator U(x) is the gauge link ordered along the infi- 
nite straight line corresponding to the "rapidity divide" 770. Both coefficient 
functions and matrix elements in Eq. (|8^) depend on the i]q but this depen- 
dence is canceled in the physical amplitude just as the scale fj, (separating 
coefficient functions and matrix elements) disappears from the final results 
for structure functions in case of usual factorization. Typically, we have the 
factors ~ (g 2 In s/m 2 — T]q) coming from the "fast" integral and the factors 
~ g 2 r]o coming from the "slow" integral so they combine in a usual log factor 
g 2 In s/m 2 . In the leading log approximation these factors sum up into the 



27 



BFKL pomeron. 

Unlike usual factorization, the expansion ( |8Cj ) does not have the additional 
meaning of perturbative versus nonperturbative separation - both the coeffi- 
cient functions and the matrix elements have perturbative and non-perturbative 
parts. This happens because the coupling constant in a scattering process is 
determined by the scale of transverse momenta. When we perform the usual 
factorization in hard (k± > fj,) and soft (fcj_ < /x) momenta, we calculate the 
coefficient functions perturbatively (because a s (k± > fj,) is small) whereas the 
matrix elements are non-perturbative. Conversely, when we factorize the am- 
plitude in rapidity, both fast and slow parts have contributions coming from 
the regions of large and small kj_. In this sense, coefficient functions and matrix 
elements enter the expansion (|80| ) on equal footing. 

3.1 High-energy OPE vs light-cone expansion 

Let me remind the idea of the usual light-cone expansion for the deep inelastic 
scattering (DIS) at moderate x. First, we take formal limit Q 2 — > oo and 
expand near the light cone (=in inverse powers of Q 2 )- The amplitude of DIS 
is then reduced to the matrix elements of the light-cone operators which are 
known as parton densities in the nucleon. At this step, the support lines for 
these operators are exactly light-like, leading to the logarithmical divergence in 
transverse momenta. The reason for this divergence is the following: when we 
expand T-product of electromagnetic currents near the light cone we assume 
that there are no hard quarks and gluons inside the proton. However, the 
matrix elements of light-cone operators contain formally unbounded integra- 
tions over fcj_, consequently there are hard quarks and gluons in these matrix 
elements. It is well known how to proceed in this case: define the renormal- 
ized light-cone operators with the integrations over the transverse momenta 
kj_ > fi 2 cut off and expand the T-product of electromagnetic currents in a set 
of these rcnormalized light-cone operators rather than in a set of the original 
unrenormalized ones (see e.g. Ref. 36). After that, the matrix elements of 

2 

these operators (parton densities) contain factors In ^ and the corresponding 

coefficient functions contain In ^ . When we calculate the amplitude we add 
these factors together, the dependence on the factorization scale [i cancels, 
and we get the usual DIS logarithmical factors ln-^. An advantage of this 
method is that the dependence of structure functions on Q 2 is determined by 
the dependence of matrix elements of the light-cone operators on fj, which is 
governed by the renormalization group. 

To get the operator expansion for high-energy scattering, we will proceed in 
the same way. At first, we take the formal Regge limit s — > oo and demonstrate 
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that the amplitude in this limit is reduced to matrix elements of the Wilson- 
line operators representing the two quarks moving with the speed of light in 
the gluon "cloud." Formally, we obtain the operators U ordered along light- 
like lines. Matrix elements of such operators contain divergent longitudinal 
integrations reflecting the fact that light-like gauge factor corresponds to a 
quark moving with speed of light (i.e., with infinite energy). The reason for 
this divergency is the same as in the case of usual light-cone expansion: the 
fast-quark propagator in the gluon "cloud" is replaced by the light-like Wilson 
line assuming that there are no fast gluons in the cloud. However, when 
we calculate the matrix element of the Wilson-line operators with light-like 
support, the integration over the rapidities of the gluon rj p is unbounded so our 
divergency comes from the fast part of the cloud which does not really belong 
there. Indeed, if the rapidity of the gluon rj p is of the order of the rapidity 
of the quark, this gluon is a fast one. As a result, it will contribute to the 
coefficient function (in front of the operator constructed from the slow fields) 
rather than to the matrix element of the operator. Similarly to the case of DIS, 
we need some regularization of the Wilson-line operator which cuts off the fast 
gluons. As demonstrated in Ref. 35, it can be done by changing the slope of 
the supporting lines. If we wish the longitudinal integration stop at i] = r] 0} 
we should order our gauge factors U along a line parallel to n — api + &P2, 
then the coefficient functions in front of Wilson-line operators (impact factors) 
will contain logarithms ~ g 2 lnl/<7. Similarly to DIS, when we calculate the 
amplitude, we add the terms ~ g 2 In 1/cr coming from the coefficient functions 
to the terms ~ g 2 In m y s coming from matrix elements so that the dependence 
on the "rapidity divide" a cancels and we get the usual high-energy factors 
g 2 In which are responsible for BFKL pomeron. Again, the advantage of 
this method is that the energy dependence of the amplitude is determined by 
the renorm-group-like evolution equations for the Wilson-line operators with 
respect to the slope of the line. 



3.2 High-energy asymptotics as a scattering from the shock-wave field. 

Consider again for simplicity the high-energy 7*7* scattering (|^). To put this 
amplitude in a form symmetric with respect the top and bottom photons, we 
make a shift of the coordinates in the currents by (z», 0, 0j_) and then reverse 
the sign of z.. This gives: 

4y e -ipA-x-ipB-y e -ia r z.+ip r z,-i(r,z) j 



A(s,t) = — i~ J d 2 z±dz m dz* j ' d A xd 

T{j A {x.,x* + z*,x± + z±)j' A (0, z*,z±) 
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x j B (y. + z.,V*,Vx)3b(z>A0±)} O). (81) 

2 

As we discussed in Sec. 1, a r ~ /3 r ~ — so it can be neglected. 

It is convenient to start with the upper part of the diagram, i.e., to study 
how fast quarks move in an external gluonic field. After that, functional inte- 
gration over the gluon fields will reproduce us the Feynman diagrams of the 
type of Fig. 3: 

A(s,t) = -t| Jd 2 z±e- i(r ' z)x Af' 1 JvAe iS ^det{iW) (82) 
'dz, I d A xe- tpA - x (Tj A (x.,x* +z*,x ± + z±)f A (0, z»,z±)) 



J dz.Jd 4 ye-^ y (Tj B (y.+z.,y*, yi _)f B (z.,0,Ox)) 



a r ' 



where 



(Tj,(x)My)) A = j v ^ ei s^A) ■ ( 83 ) 

Here S(A) and S(ip, A) are the gluon and quark-gluon parts of the QCD action 
respectively, and det(iV) is the determinant of Dirac operator in the external 
gluon field. 

The Regge limit s — > oo with p A and p B fixed corresponds to the following 
rescaling of the virtual photon momentum: 

PA = A^ 0) + — J§ pa, (84) 

2\ P r ■ p 2 

with ps fixed. This is equivalent to 

Pl = Ap<°\ P 2=P 2 °\ (85) 

where p^ and p^ are fixed light-like vectors so that A is a large parameter 
associated with the center-of-mass energy (s = 2Xp^' -p^)- Let us study the 
asymptotics of high-energy 7*7* scattering from the fixed external field 

dx J dz6(z.)e-^ x -^ z ^ (T{j^(x + z)j v {z)}) A . (86) 

Instead of rescaling of the incoming photon's momentum (]84|), it is convenient 
to boost the external field instead: 

dxdz5(z.)e-* PAX -^ z ^ (T{j^x + z)j v {z)}) A 
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dxdz6(z )e- i ^ *-« r >^ (T{j^x + z)j v (z)}) B , (87) 

where p^ 1 = p^ + ^P2 and the boosted field has the form 

B (x ,x*,x±) = A^4o(^,a;*A,a;j_), 
Btf (x , Xtf , x±) = — A* (-^— , A, x± ), 
Sx(x ,a;*,a;x) = Ai(j, sc*A, (88) 

where we used the notations x Q = x^p^, x* = x^pi^. The field 

Au(x , iejJ = XoPj - 1 H + sex) (89) 

so s 

is the original external field in the coordinates independent of A, therefore we 
may assume that the scales of x ,x* (and x±_) in the function (^9|) are O(l). 
First, it is easy to see that at large A the field B fJl (x) does not depend on x a . 
Moreover, in the limit of very large A the field B^ has a form of the shock 
wave. It is especially clear if one writes down the field strength tensor G^u 
for the boosted field. If we assume that the field strength for the external 
field vanishes at the infinity we get 

Goi(x , x*, x±) = \F i(-?-,x*\,x±) — > 5(x*)Gi(x±), 
1 x 

G ri (x ,x*,xj_) = -F ti (-^,x*\,x±) — * 0, 

G *(x , x„, x±) — F *(-^-,x*\,xx) — * 0, 

A 

G ik (x ,x*,x±) = F ik (-j-,x*\,x±) — > 0, (90) 

so the only component which survives the infinite boost is -F x and it exists 
only within the thin "wall" near = 0. In the rest of the space the field B^ 
is a pure gauge. Let us denote by Q the corresponding gauge matrix and by 
B n the rotated gauge field which vanishes everywhere except the thin wall: 

B? = lim ||-G? o (0, \x*,x x ) - 6(x*)^G?(x x ), = B ± = 0. (91) 

To illustrate the method, consider at first the propagator of the scalar 
particle (say, the Faddeev-Popov ghost) in the shock-wave background. In 
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Schwinger's notations we write down formally the propagator in the external 
gluon field A fl (x) as 



G(», V ) = ((a 



p2 



(p + + it 



(92) 



where ((x|y)) = S^ 4 \x — y), 

((^Kb)) - -i-^5W(x-v), ((x\A,\y))=A,(x)5^(x-y). (93) 

Here |a;)) are the eigenstates of the coordinate operator X\x)) = x\x)) (nor- 
malized according to the second line in the above equation) . From Eq. ( p3| ) it 
is also easy to see that the eigenstates of the free momentum operator p are 

\x)). The path-integral representation of a 



-ip-x 



the plane waves \p)) = J d 4 
Green function of scalar particle in the external field has the form: 

1 



V 2 

OO 



y =-* 



dr 


x(r)— x 



((• 



JtV 



(94) 



Vx(t)e~ l fo dt -pexp{ig I dt{B™ (x(t))x» (i)} 



where r is Schwinger's proper time. It is clear that all the interaction with 
the external field occurs at the point of the intersection of the path of 
the particle with the shock wave (see Fig. [13). Therefore, it is convenient to 




Figure 13: Propagator in the shock-wave field. 
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rewrite at first the bare propagator 



((-> 



47r 2 (x — y) 2 



x(r)—x 



x(0)=y 



Vx(t)(r)e~ 



Jo 



(95) 



marking the point of the intersection of integration path with the plane = 0. 
To this end, consider the case x* > 0, < and insert 



1= (It'x*(t')5{x*{t')) 



(96) 



in the path integral (95). (Here r' has the meaning of the time at which the 
intersection with the plane z* = takes place). We get 



1 



i I driix 



(97) 



rfr / dr'N- 1 
ii Jo 



x(t)=x 



(o)=y 



C Vx(t)x*{T')6{x*{T'))e-' 1 fo dti 



= -i dr dr' I dzSiz^Af- 1 
Jo Jo J 

X(T') = Z r. T > 2 

Vx(t)x,{T')e~ l Jo dt ~ 



x(t)=x 



Vx 



(i>- ,; J>^Ar 



x(0)=y 

Making the shift of integration variable r — t' — > r, we can rewrite the path 
integral (|97j) in the form: 



/>oo />oo p 

—if dr dr 1 / dz5(z st ) 
Jo Jo J 

Af- 1 vx(ty l h dt ^N- x \ 

Jx(0)=z Jx 



(98) 



x(r') = 



c(0) 

Using Eq. ( |95| ) and similar formula 

r x(r) 



(0)=y 



Vx(t)x*e 



. -i r dti 



dTjV~ 



(0)=y 



Vx{t)x^{T)e Jo * = — £ 



7r 2 (x — y) 4 ' 



(99) 



we arrive at the following representation of the bare propagator (in the case of 
x* > 0, y* < 0): 



1 



jr + ie 



dzS(z* 



1 



47r 2 (x-z) 2 7r 2 (z-t/) 4 



(100) 
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where z is the point of the intersection of the path of the particle with the 
shock wave. 

Now let us recall that our particle moves in the shock-wave external field 
and therefore each path in the functional integral is weighted with the 

additional gauge factor Pe g J B ^ dx ^ . Since the external field exists only within 
the infinitely thin wall at = we can replace the gauge factor along the 
actual path x^(t) by the gauge factor along the straight- line path shown in 
Fig. It intersects the plane = at the same point (z Q ,z±) at which 
the original path does. Since the shock-wave field outside the wall vanishes 
we may formally extend the limits of this segment to infinity and write the 
corresponding gauge factor as U (z±) = [oop% + z±, — oopi + z±\. The error 
brought by replacement of the original path inside the wall by the segment of 

straight line parallel to pi is \J~^-- Indeed, the time of the transition of the 
particle through the wall is proportional to the thickness of the wall which is 

2 

~ — . It indicates that the particle can deviate in the perpendicular directions 
inside the wall only to the distances J Thus, if the particle intersects this 

wall at some point (z*,z±) the gauge factor Pe 9 J B ^ dx " reduces to U n (z±). 
One can now repeat for the path integral (194) the steps which lead us from 



path-integral representation of bare propagator ( |95| ) to the formula (100); the 
only difference will be the factor U (z±) in the point of the intersection of the 
path with the plane = 0: 



1 



dz6(z*) 



1 



4ir 2 (x — z)' 



:U Q (z X ) 



V* 



(101) 



(in the region x* > 0,y* < 0). It is easy to see that the propagator in the 
region < 0, y* > differs from Eq. (101) by the replacement U n ^ U n '. 
Also, the propagator outside the shock- wave wall (at < or > 0) 

coincides with the bare propagator. The final answer for the Green function 
of the scalar particle in the B Q background can be written down as: 



V- 



(x*j/*) + / dz8{z*) 



47r 2 (x-y) 2 v " ' J y '4ir 2 {x-z) 2 
x {U n (z ± )e(x*)6(-y*) - U^(z ± )6(y«)6(-x«)} 



y* 



(102) 



^(z-yY 



We see that the propagator in the shock-wave background is a convolution of 
the free propagation up to the plane = 0, instantaneous interaction with 
the shock wave described by the Wilson-line operator U n (U^), and another 
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free propagation from z to the final point (see Fig. |l^) One can check that the 
Green function ( |102| ) is continuous as — * (or — ► 0). 

In order to get the propagator in the original field we must perform 
back the gauge rotation with the Q, matrix. It is convenient to represent the 
result in the following form: 



i x \h\ y l 



where 



[x,y]0(x*2/») + / dz6(z») 



1 



4Tr 2 (x-y) 2L ' aJ v a ' J v '4tt 2 (x-z) 2 
x{U(z x \x,y)0{x.)O(-y m ) - rt(z±;x,y)0(y m )9(-x.)} 



(103) 



tj* 



ir 2 (z- y y 



U{z±;x,y) = [x,z x ][z x ,Zy][z y ,y] , 
2 s 

— X*p 2 ,Zj_), Zy 



I 2 (o) 
{—z p\ 

so 



•S'C) 



y* 



(104) 



is a gauge factor for the contour made from segments of straight lines as shown 
in Fig. |l4[ Since the field outside the shock-wave wall is a pure gauge, the 
precise form of the contour does not matter as long as it starts at the point x, 
intersects the wall at the point z in the direction collinear to p2, and ends at 
t he p oint y. We have chosen this contour in such a way that the gauge factor 
( |l04| ) is the same for the field B^ and for the original field (see Eq. (f^sj)). 

The quark propagator in a shock-wave background can be calculated in a 
similar way (see Appendix 7.2), 



2tt 2 (x - y)' 



[x, y]6{x+y+) + i dz5(z*) 



2tt 2 (x-z) 4 



x {U(z±;x, y)9(x* )8(~y* )-U\z ± ;x, y)9(y* )fl(-s«)} 2 J { ~^ y)4 ■ (105) 
For the quark-antiquark amplitude in the shock-wave field (see Fig. [u|) we get 



(106) 



1 






l 


r 


y)i»{ 


h 


T 



Tr 7 ,(( 

Tr7 M (^ - j)^ v {j - f) 



x Tr 7p 



47r 4 (x - y) 8 



0(x*y* 



2tt 2 (a; - z) 4 ^ 2 2tt 2 (z - y) 4 ' 2tt 2 (y - z') 4 ^ 2tt 2 {z' - x) 4 



Q{— x*y*) / dzdz' 8{z*)5{z' tl 



i> 2 



U(zx;zi), 



where we can write down the gauge factor U(z±; z' ± ) = U(z±; x, y)W(z' ± ; y, x) 
as a product of two infinite Wilson- lines operators connected by gauge segments 
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Figure 14: Quark-antiquark propagation in the shock wave. 



at ±00, 

= lim < [upi + zj_,—upi + z±][— upi + z_l, —upi + z' ± ][— upi + z'_ L ,up\ + z'A 
x[upi + z' x ,up 1 + z ± ]j = U z [z±,z' x }-Ul,[z' ± ,z_ L ] + , (107) 

Here we use the notations 



-oopi + zj_, — oopi 



[x±,yx]+ = [oopi +zj_,oopi +z' ± ], [x±,y±\- 

. (1 ° 8) 

As we mentioned above, the precise form of the connecting contour at infinity 
does not matter as long as it is outside the shock wave. We have chosen this 
contour in such a way that the gauge factor (107) is the same for the field 
and for the original field (see Eq. (p8[)). Now, substituting our result 
for quark-antiquark propagation ( 106 ) in the right-hand side of Eq. (|8rj), one 
obtains 



d 4 x / d 4 z6(z.)e- l ^ z ^e-^ x (T{j A (x + z)j' A {z)}) A 



d 2 k 



-L j A 



4tt 2 



I A (k ± ,r ± )Ti{U(k ± )U^r ± -k ± )}, (109) 
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where the impact factor I A is given by Eq. (15). For brevity, we omit the end 



gauge factors (108) 



Formula (107) describes a quark and antiquark moving fast through an 
external gluon field. After integrating over gluon fields in the functional inte- 
gral we obtain the virtual photon scattering amplitude (p2|). It is convenient 
to rewrite it in the factorized form: 

A( PA ,p B )^ l S -Y / ^ J ^I A (k ± ,r ± )((Tr{U(k ± )U\r ± -k ± )})). (110) 

where I A (p±) = e A e A I A (px)- The gluon fields in U and W have been pro- 
moted to operators, a fact which we signal by replacing U by U, etc. The 
reduced matrix elements of the operator Tr{U(kx)U' (r± — k±)} between the 
"virtual photon states" are defined as follows: 

((Tv{U(k ± )U\r ± - k ± )})) = Jd 2 x ± e-^H(^{U(x±)UH0)})) 

((Tr{U(x ± )U\0)})) = - Jd 4 z5(z«y ir < z)± Jd 4 y e - tpB y (111) 

(0\T{Tr{U(x ± )tf(0)}j B (y + z)j' B (z)}\0). 

This matrix element describes the propagation of the "color dipole" in the 
background of the shock wave created by the second virtual photon. 

It is worth noting that for a real photon our definition of the reduced 
matrix element can be rewritten as 

(e,p B \Tr{U(x ± )tf(x' x )}\e',p B +(3p B ) = 2ttS((3) ((Tr{U(x ± )& (x' x )})), 

(112) 

where e and e' represent the polarizations of the photon states. The fac- 
tor 2-k8{0) reflects the fact that the forward matrix element of the operator 
U(x±)W (x' ± ) contains an unrestricted integration along p\. Taking the inte- 



gral over (3 one reobtains Eq. (Ill) 



3.3 Regularized Wilson-line operators 



In the Regge limit (84) we have formally obtained the operators U ordered 
along the light-like lines. Matrix elements of such operators contain divergent 
longitudinal integrations which reflect the fact that light-like gauge factor cor- 
responds to a quark moving with speed of light (i.e., with infinite energy). This 
divergency can be already seen at the one-loop level if one calculates the con- 
tribution to the matrix element of the two- Wilson- line operator U(x±)W (y±) 
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between the "virtual photon states" . As I mentioned above, the reason for this 
divergence is that we have replaced the fast-quark propagators in the "external 
field" represented by two gluons coming from the bottom part of the diagram 
in Fig. [l5| a by the light-like Wilson lines in Fig. [T^b. The integration over 




(a) (b) 



Figure 15: A typical Feynman diagram for the 7*7* scattering amplitude (a) and the corre- 
sponding two- Wilson-line operator (b). 

rapidities of the gluon rj p in the matrix element of the light-like Wilson-line 
operator U(x±)U' (y±) is formally unbounded , consequently we need some 
regularization of the Wilson-line operator which cuts off the fast gluons. As 
demonstrated in Ref. 35, it can be done by changing the slope of the supporting 
line. If we wish the longitudinal integration stop at r\ = 770, we should order 
our gauge factors U along a line parallel to Pq — pi + Cpz where 770 = ln£. e 
We define 

U^(xj_) = [dOp^ + Xj_, — OOp C + 5Ej_], 

e The situation here is again quite similar to the usual OPE for DIS. Recall that when 
separating the Feynman integrals over loop momenta p into the coefficient functions (with 
p 2 2> fJ?) and matrix elements (p 2 -C £t 2 ) we expand hard propagators in powers of soft 
external fields. As a result of this expansion we formally obtain the expressions of the type 
V>(Aei)[Aei, 0]i/>(0) with external fields lying exactly on the light cone. In operator language 
it corresponds to the matrix element of the same light-cone operator i/>(Aei)[Aei , 0]?/>(0) 
normalized at the point fi 2 in order to ensure the restriction that matrix elements of this 
operator do not contain virtualities larger than /x 2 . Moreover, in principle we can regularize 

these light-cone operators for DIS by changing the slope of the supporting line (say, take 
2 

e = ei + 7^62). The only reason why we use the regularization by counterterms is that, 
unlike the regularization by the slope, counterterms are governed by renormalization-group 
equations. 
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U'^(x±) = [—oop^-\-xx,oop^ + xx\- (113) 
Matrix elements of these operators coincide with matrix elements of the oper- 
ators U and ffl calculated with the restriction a < a = \/^4 imposed in the 

V s s 

internal loops (and external tails). Let us demonstrate this using the simple 
example of the matrix element of the operator U**(kx)U' [rx — kx) coming 
from the diagram shown in Fig. [l5]. It has the form 

* 6 fda p d 4 p' [jap - 2c/ k )p k s - {k + k')\]^ B {k') 



2 " 3 -I L>n iiy " 4 (Ca2 s + % - ie) 2 (a' fc /3^ - + ie) 2 

' (H4) 



-{ctp - a')(a p C + (3' k )s - (k - k') 2 x + i 



where the numerator comes from the product of two three-gluon vertices (18) 

4r*. CT (fc,-fc')r*. ff (fc, -k') == K - 2a' fe Ks - {k + k')\. (115) 

As we shall see below, the logarithmic contribution comes from the region 
v , > a k > o/ fe - 1 > /3' p > p = -(a k ~ \J~^- In this region 
one can perform the integration over (3' k by taking the residue at the pole 

-(a p - a')(a p ( + j3' k )s - (k - k')\ + it . The result is Q 
- 6 r da k da'u f d 2 k\ 



•s 



-± [e(a k > a' k > + 6(0 > a' k > a k )] (116) 



2?r 2vr J Air 2 
(fc 2 + p'l al( S /2) $ s (a' kPl (a k ( + %|^)p 2 + k' ± 

\a p - a' k \ (Cap + k \ - it) 2 [^{k-k')l +p'\ + ze] 2 

We see that the integral over a p is logarithmic in the region J ^ a p 

a' k ~ — (cf. Eq. (18)). The lower limit of this logarithmical integration is 
provided by the matrix element itself ((3 k ~ 1 in the lower quark bulb) while 
the upper limit, at a 2 <~ m 2 /Qs is enforced by the non-zero C an d the result 
has the form 

((TrUHkx)U^(-kx))) Fis . 15 = fyn(-^) [ f^%±^I^k' x ). 

8tt \m 2 CJ J 4tt 2 fc 4 p ,* 

(117) 



•'In the region we are investigating, we can neglect the /3'. dependence of the lower quark 
loop. 
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Similarly to the case of usual light-cone expansion, we expand the ampli- 
tude in a set of "regularized" Wilson- line operators tK (see Fig. pj| ): 

A(p A ,PB =>■ p'aiPb) = ^ J d 2 x 1 ...d 2 x n C(x 1 , ...x n : () (118) 
x {p B \TT{U^x 1 )U K {x 2 )...U^{x n ^)U^{x n )\p' B ). 
The coefficient functions in front of Wilson-line operators (impact factors) will 





+ 



(a) 




Figure 16: Decomposition into product of coefficient function and matrix element of the 
two- Wilson-line operator for a typical Feynman diagram. (Double Wilson line corresponds 
to the fast-moving gluon.) 

contain logarithms ~ g 2 In 1 J a and the matrix elements ~ g 2 In . Similar to 
DIS, when we calculate the amplitude, we add the terms ~ g 2 \n\/a coming 



from the coefficient functions (see Fig. |16|b) to the terms ~ g In m i/ s coming 
from matrix elements (see Fig. |l6|a) so that the dependence on the "rapidity 
divide" a cancels resulting in the usual high-energy factors g 2 In -Ag which are 
responsible for the BFKL pomeron, cf. (|50]). 



In the LLA, the light-like operators U and W in Eq. (110) should be re- 
placed by the Wilson-line operators (ft and (ft^ ordered along n \\ pa- Indeed, 
let us compare the matrix element (117) shown in Fig. 6b to the corresponding 
physical amplitude ( |l7|) shown in Fig. 6a. The int egral in Eq. ( |l7| ) is similar 
to the one for the matrix element of the operator ( 117 ), except that there is 
now a factor of the upper quark bulb and the integral over p± . If we calculate 
only the contribution of the diagram in Fig. 6a , we would get (cf. Eq. (pq)) 



47T 



d 2 k ± d 2 k' ± k\ + p> ± 



47T 2 47T 2 



kip'i 



I A (k ± )I B (k' ± ) (119) 



which agrees with the with estimate Eq. ( |117| ), if we set ( — This cor- 
responds to making the line in the path-ordered exponential collinear to the 
momentum of the photon. 
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3.4 One-loop evolution of Wilson-line operators. 



As we demonstrated in previous section, with the LLA accuracy, the improved 



version of the factorization formula Eq. (109) has the operators U and W 

2 

"regularized" at Q ~ — : 



d 4 x I d 4 z5(z.)e- lpA - x -^ r ' z)± T{j A (x + z)j' A (z)} 
d 2 k A 



(120) 



«L rA 



4tt 2 



I A {k±, rj.)Tr{[/ c =^ (k)U K=n ^ (r - k)}. 



In the next-to-leading order in a s we will have the corrections 
~ a s TrU(x x )UHy±) r TrU(y ± )UHz 1 _), see Fig. [l§ 

Next we derive the equation for the evolution of these operators with 
respect to slope £ (in the LLA). In order to find the behavior of the matrix 
elements of the operators (x±)U'^ (y±) on the slope C we must take the 
matrix element of this operator "normalized" at (4 and integrate over the 



momenta with o\ 



> a > 02 



(similar to the case of ordinary 



Wilson OPE where in order to find the dependence of the light-cone operator 
on the normalization point [L we integrate over the momenta with virtualities 
fif > p 2 > fJ%). The result will be the operators U and W "normalized" at the 
slope (2 times the coefficient functions determining the kernel of the evolution 
equation. The calculation of the kernel is essentially identical to the calculation 
of the impact factor with the only difference of having initial gluons instead of 
quarks. Here we will present only the outline of the calculations; the details 
can be found in Appendix C. 

In the first order in a s there are two one-loop diagrams for the matrix 
element of operator [/(xjjf/t [yj_) in external field (see Fig. [It]). This external 





.-• ., r 



Figure 17: One-loop diagrams for the evolution of the two- Wilson-line operator. 

field is made from slow gluons with a < (2- Like the case of the fast quark 
propagator considered above, it is convenient to go to the rest frame of "fast" 
gluons, as a consequence the "slow" gluons will form a thin pancake. 



41 



Let us start with the diagram shown in Fig. p"7|a. We will calculate the one- 
loop evolution of the operator U(x±) <g>ffi (y±) = {U(x±)}j{W(y±)}f with the 
non-convoluted color indices. In the LLA, the slope of the operators U can 
be replaced by pi . Using the expression for the axial-gauge gluon propagator 
in the external field (13081) f\ we obtain: 



(121) 



ig 2 / du[oopi,upi] x t a [upi,~copi] x / dv[-oopi, vpi] y t b [vpi, oopi] y 



x ((u Pl + x± ( Plt - V.^)0^(p lv - -Hi-V.) 
W p ■ pi P-P2 

Hereafter we use the space-saving notation 

[un, vn] x = [un + xj_, vn + xj_] 



vpi + y± 



ab 



(122) 



We may drop the terms proportional to V % in the parenthesis since they 
lead to the terms proportional to the integrals of total derivatives, namely 

J du[oopi,upi}t a [upi,-oopi]pi tJ/ {D^^(upi, ...))ab 

/d 
du— {[oopi,upi]f a [wpi,-oopi]($(Mpi, ...))ab} = (123) 

and similar for the total derivative with respect to v. Now, we can rewrite 
Eq. ([i2ll) in the form 
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du[oop 1 ,up 1 ] x t a [upi,-oop 1 ] x (124) 



dv[-copi,vpi]yt b [vpi,copi}y(^upi +x± O,, vpi +V±Jj ■ 

As in the calculation of the quark propagator, it is convenient to go to the 
rest frame of "fast" gluons. In this frame the "slow" gluons will form a thin 
pancake shown in Fig. |l8|. At first, we consider the case > 0, y* < 0. It is 
clear from the picture that we can rewrite Eq. (|l24|) as follows: 



ig 2 t a U(x ± )®t b U\y ± ) 



(125) 



du I dv((up^ 



■Xj_ 



a. 



(0) 

vp A 



y± 



a b 



,t, fii rthnr terms in expansion in powers of gluon propagator ( 306 ) 



9 It can be demonstrated th. 
beyond those given in Eq. (307) do not contribute in the LLA. 
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a b c 

Figure 18: Path integrals describing one-loop diagrams for Wilson-line operators in the 
shock-wave field background. 



(we shall calculate only the contribution ~ U which comes from the region 
> 0, < - the term ~ W coming from x* > 0, y* < is similar). 
Technically it is convenient to find at first the derivative of the integral of 
gluon propagator in the right-hand side of Eq. ( |125| ) with respect to x± . Using 
the thin-wall approximation we obtain 



a. 



s 2 f, x( Anjx- zf 
2 J 167r z x* 

x {2[FF](z ± )-i[DF](z ± )} 



(126) 



1 



4tt 2 (z - yY 



where 



[DF](x ± ) 
[FF](x x ) 



dcf 



del' 



du[oopi, upi] x D a F a .(up 1 + x±)[upi, -oopx] x , 

du J dvQ{u — v)[copx, up\] x F^(upx + x±) 
[up 1 iVpi\ x F(..(vpi + x±)[vp!,-oopi] x . (127) 



It is easy to see that the operators in braces are in fact the total derivatives of 
U and U' with respect to translations in the perpendicular directions, 



d 2 ± U\x ± ) 



dxidx. 

d 2 
dxi dx 



■U(x x ) = -i[DF](x x ) + 2[FF](x ± ), 
■U\xa_) = i[DF](x±) + 2{FF][x±), (128) 



(note that djU = -d 2 U). 
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For the derivative of the gluon propagator (x\piO\y) we obtain: 



-w 



du / dv 



16tt 4 



dz i 



(o) , 

UP A + 

00 du , 
— dv 



dz. 



(0) 



2/-L 



(129) 



(.t_l - z ± _) l [d\U{z 1 _)] ab 



[u(u(s — 2z.) - (x - z)\ — ie][v{vQs + 2z.) — (y — z)\ — ie] ' 
The integration over z, can be performed by taking the residue; the result is 
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16tt 3 



dz± I — dvj-—r 

u [{x - 



(x ± - zx) l [dj_U(zx)] a b 



z)j_ v (v ~ z)±_ v ~ uv(u + v)( > s + ie] 



(130) 



1 x ' vy "jx l 

This integral diverges logarithmically when u — > — in other words when the 
emission of quantum gluon occurs in the vicinity of the shock wave. (Note 
that if we had done integration by parts, the divergence would be at v — > 0, 
therefore there is no asymmetry between u and v). The size of the shock wave 
Zst r~j m^ 1 ^ (where 1/m is the characteristic transverse size) serves as the 
lower cutoff for this integration and we obtain 
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16tt 3 

„2 



0-2 



dz± 



1 da (x_l- z_ L ) i [dj_U(z ± )] ab 
a [{x — z)\a + {y — z) 2 ^] 



r In 

167I" 3 (72 



(131) 



(recall that a = 1 — a). Thus, the contribution of the diagram in Fig. 
the LLA takes the form 




(U(x ± )W(y ± )) A = - 

+ U(x ± )t a ®U^(y±)t b 



t a U(x ± )®t b U\y ± ] 



Xj_ 



1 

p± 



Px 



y± 



i(d 2 xU^ 
Px Px 



y± 



(132) 



where we have added the term coming from x* < 0, > 0. A corresponding 
result for the diagram shown in Fig. |l8| b can be obtained by comparing the 
space-time picture Fig. 18 3 for this process with Fig. Km,, 



{Ut(xx)®UK(yx)), 



^In^l 

27T (72 



U(xx)®t a U\yx)t b 



y± 



Px Px 



y± 



(133) 
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Likewise, the diagram in Fig. [18p yields 



((- 



p± p± 



(134) 



The total result for the one-loop evolution of two- Wilson-line operator is 
the sum of Eqs. jl32|), (pf), and (pi, 



({^(^)}}{c> tCl (y±)}f)A = ^u^Jdz ± 



(135) 



(x-z,y- z)j_ 
(x-z) 2 x (y- z)\ 



{U^(z ± )})Tr{U^(x ± )U^(z ± )} - N c {U^{x±)}<\ {U K2 (y±)f } 

{u^(x ± )}\ \{u K Hz±)}i^{u C2 (z±)u K2 (y±)} 



(x-z) 2 x 



N c {U K2 (y ± )} 



The evolution of a general n- Wilson-line operator is presented in Appendix 
7.3.Q 

3.5 BFKL pomeron from the evolution of the Wilson-line operators 

As we demonstrated in Sec. 3.2, with the LLA accuracy the improved version of 
the factorization formula Eq. (109) has the operators U and W "regularized" 

atC~4: 



otx I d 4 zS(z.)e- lpA - x ' l ^ z ^T{j A (x + z)f A (z)} 
d 2 k. 



(136) 



E 



4tt 2 



/ A (fc_L, r ± )Tr{[/ c =^ {k)U K =^ (r - k)} + 0(g 2 ). 



h A more careful analysis performed in Appendix shows that the Wilson lines U and U< 
connected by gauge links at infinity, see Eq. (E99). 
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In the next-to-leading order in a s we will have the corrections 

a s TrU(x ± )W (y ± )TrU(y ± )W (z±), see Fig. pi The matrix element of this 



operator ((U^(x±)U^(y±))) (see Eq. (Ill) for the definition) describes the 



gluon-photon scattering at large energies ~ s. (Hereafter we will wipe the 
label (") from the notation of the operators). The behavior of this matrix 
element with energy is determined by the dependence on the "normalization 
poi nt" C - From the one-loop results for the evolution of the oppr^tors U and 
W (|135|) it is easy to obtain the following evolution equationc 



17/ 



(-Q^U(x±,y±) = -^^■Jdz±^U(xx,zx)+U(zx,y±)-U(xx,y±) 

(g ~ v) 2 x 
(xx - zx) 2 {zx - iJx) 2 ' 



U(x,z)U(z,y)\- { *j} 2± _, 2 , (137) 

J (x i - z \ Viz i v \ ) 



where 



U(xx,yx) = ^-(Tr{U(xx)[xx,yx}-UHyx)[yx,xx]+} - N c ) (138) 



(cf. Eq. (107)). Note that right-hand side of this equation is both infrared 
(IR) and ultraviolet (UV) finite^] Wc see that as a result of the evolution, 
the two- line operator Tr{UU>} is the same operator (times the kernel) plus 
the four-line operator Tr{UW}Tr{UW}. The result of the evolution of the 
four-line operator will be the same operator times some kernel plus the six-line 
operator of the type Tr{UW}Tr{UW}Tr{UW} + Tr{UWuW}Tr{UW} and 
so on. Therefore it is instructive to consider at first the linearization of the 



Eq. (137) with the number of operators U conserved during the evolution. 



The linear evolution of the two- line operator U(xx,yx) is governed by the 



'The similar non-linear equation describing the multiplication of pomerons was suggested in 
Ref. 38 and proved in Ref. 39 in the double-log approximation 

J The IR finiteness is due to the fact that TrC/f/^ corresponds to the colorless state in t- 
channel, as a consequence the IR divergent parts coming from the diagrams in Figs. |l8|a, 
Ob, and Est cancel out. Tf we had the exchange by color state in t-channel, the result will 
be IR divergent (cf. Eq. (uq)). 
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BFKL equation [] 

Q^-U{xx,Vx) (140) 

= ~4^5 N ° J dz± {^(^-L^-l) + U( z x,yx) -U(x±,yx)} (g_ (JZ ■ 

Let us start from the simplest case of forward matrix elements (which descr ibes , 
for example, the small-x DIS from the virtual photon). Then the equation (140) 
takes the form 



C^«W(*l)» = -^N c Jdz ± [U(x - z ± ) +U(z ± )-U(x ± )] X ^ 



(141) 



where ((U(x±))) = ((U(x±,0))) (see Eq. ([111])). The eigenf unctions of this 
equation are powers (xj_)~^ ,v and the eigenvalues are —^N c x{v), where 
x(u) = — Re^(| + iv) — C. Therefore, the evolution of the operator U takes 
the form: 



x fdz ± (^)-i-^'(z))) (142) 



We may proceed with this evolution as long as the upper limit of our 

logarithmic integrals over a, \J ^A, is much larger than the lower limit ^f- de- 
termined by the lower quark bulb, see the discussion in Sec. 3.3. It is convenient 
to stop evolution at a certain point £o such as 

(o = u 2 ^, a«l, .g 2 l n( 7«l, (143) 

then the relative energy between the Wilson-line operator and lower vir- 
tual photon will be sq = m 2 a 2 which is big enough to apply our usual high- 
energy approximations (such as pure gluon exchange and substitution g^ v — > 

fe If F(k±,r±) satisfies the BFKL Eq. @ then 

U ix±>y±) = /^e^.W-^ (139) 



( ™ k± ' r ±> - l[S(k±) + S(rx - kx)} dk' 
\k\{f-k)\ 1 J ± (k>)l(r-k>) 
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j^P2fiPiv) but small in a sense that one does not need take into account the 
difference between g 2 In and g 2 In m 2 a2 ■ Finally, the evolution ( |140D takes 
the form: 



J 2ir z m z J 

, (144) 

Now let us rewrite this evolution in terms of original operators l/W in the 
momentum representation. One obtains: 



((Tt{U^( P± )U^^(- P± )})) 



dv 
2k 2 



(f±y 



(m 2 ) 



(145) 



where we omit the gauge links at infinity ( 108 ) for brevity. Since we neglect the 
logarithmic corrections ~ g 2 In a the matrix element of our operator U^ U^ a 
coincides with impact factor I B up to 0(g 2 ) corrections: 



(<TV{^°(px)£/ tCo (-f>±)}>) 



(146) 



*N 2 



da $ B (a p p! - ( a p p 2 + p±) 



(Coa 2 +0i) 2 



Combining Eqs. (jHOj) , (145), and (|l4q ) we reproduce the leading logarithmic 
result for virtual 77 scattering (|59|). 

In the case of small-x DIS from the nucleon the matrix^ element of the 
operator UU* describes the propagation of the "color dipole'EO in the nucleon. 
The evolution of the matrix element (N\U\N) is the same as Eq. (145) with the 
only difference that the lower impact factor I B should be substituted by the 
nucleon impact factor I N determined by the matrix element of the operator 
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UW between the nucleon states :| 
(N,p B \Tr{Ut°(x ± )UK°(0)}\N,p B +j3p 2 ) = 2n5(f3) f ^± e ^^I N { Pl _) 



Air 2 p\ 

(147) 

where 27r5(/3) reflects the fact that matrix element of the operator UW con- 
tains unrestricted integration along p^° , (cf. Eq. (111)). The nucleon impact 



factor I B (p±) defined in (147) is a phenomenological low-energy characteristic 
of the nucleon. In the BFKL evolution it plays a role similar to that of a 
nucleon structure function at low normalization point for DGLAP evolution. 
In principle, it can be estimated using QCD sum rules or phenomenological 
models of nucleon. 

In conclusion, let us present the results for the linear evolution for the 
non-forward case. Due to the conformal invariance of the tree-level QCD the 
eigenf unctions of the equation (|140|) are powersE3 



(148) 



(x-y)\ 



(x - xo) 2 ± (y- x ) 2 ± 



where xo is arbitrary. The eigenvalues are the same as for the forward case, 
— ^-N c x{v). The corresponding formula for the result of the evolution of the 
two- Wilson-line operator has the form: 



U^{x^ Vl _) = 



/ dvd *v [<rmrm 

f A \ — ^JVcX(f)) 



U C2 (x ,v) (149) 



where 



(150) 



'This is called "hard pomeron" contribution to the structure functions of DIS since the 
transverse momenta in our loop integrals are large (~ Q 2 ), at least in the lowest orders in 
perturbation theory. However, due to the diffusion in transverse momenta the characteristic 

size of the in the middle of gluon ladder is Q 2 e~ v's 2 ln s (see the discussion in Sec. 2.4), 
so at very small x the region p± ~ ^QCD ma y become important. It corresponds to the 
contribution of the "soft" pomeron which is constructed from non-perturbative gluons in our 
language and must be added to the hard-pomeron result. 
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It is worth noting that at large momentum transfers — t = fj_ m 2 N the 
nucleon impact factor is determined by the well-studied electric and magnetic 
form factors of the nucleon 

I N (k u r±) P± = m2 6 xx ,F? +n (t) + l^-uip'A') A fMP,m +n (t), (151) 

which gives an opportunity to calculate the amplitude of deeply virtual Conm- 
ton scattering from the nucleon at small x without any model assumptions e!j 



3.6 Non-linear evolution of Wilson lines 

Unlike the linear evolution, the general picture is very complicated: not only 
the number of operators U and IP increase after each evolution but they form 
increasingly complicated structures like those displayed in Eq. (153) below. In 
the leading log approximation the evolution of the 2n-line operators such as 
Tr{UW}Tr{UU'}...Tr{UW} comes from either self-interaction diagrams or 
from the pair-interactions ones (see Fig. FL9J) The one-loop evolution equations 



::: ::: ::. 




r 



• ? : r 

< »•* < »•? 



B 



Figure 19: Typical diagrams for the one-loop evolution of the n-line operator. 



for these operators can be constructed using the pair-wise kernels calculated in 
the Appendix C. For instance, the evolut ion e quation for the four-line operator 
appearing in the right-hand side of Eq. ( 137 ) has the form: 



C^-Tr{U x [x, z]-Ul[z, x] + }Tt{U z [z, y]-U$\y, z} + } 



(152) 



16?r 3 



dt\ 



Tr{U x [x, t\-Ul [t, x} + }Tr{U t [t, z]-Ul\z, t]+} 
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N c Tr{U x [x,z]-U$[z,x} + }]Tr{U z [z,y]-Ut{y,z] + }- U ' ' 



(x-t)l(z-?)i 



+ Tr{U x [x,z].Ut[z,x]+}- {V Z)2± 



x Tr{U z [z, t].U}[t, z}+}TY{U t [t, y]-U}\y, t}+} 

- N c Tr{U z [z,y]-Ul[y,z}+} 

+ [Tr{U x [x, z]-Ut[z, t]+U t [t, y]-Ul[y, z]+U z [z, y]-U}[t, x}+} 
+ Tr{U x [x, t] _ U} [t,z} + U z [z,y\- U$ [y,t} + U t M - U\ [z,x]+} 

- 2Tr{U x [x,y]-Ul}[y,x]+ 

1 



(x-t,y- f)_L 
(x-?) 2 x (y-t) 2 x (z-t) 2 

(S-t,z-t)x {z-t,y-t)± 



(x-t) 2 ± (z-t)l (z-t)i(y-?)l 



where we have displayed the end gauge links ( f299| ) explicitly. Note that each 
of the separate contributions (30C) and (301) corresponding to the diagrams in 
Fig. |39|a and |39]b diverges at large t while the total result (152) is convergent. 



This is the usual cancellation of the IR divergent contributions between the 
emission of the real (Fig. |39|a) and virtual (Fig . |39"fc) gluons from the colorless 
object (corresponding to the l.h.s. of Eq. (152)) (cf Eq. ( 137 )). 

Thus , the result of the evolution of the operator in the right-hand side of 
Eq. (12C) has a generic form: 

Tr{Ut[x,y}-UK[y,x} + }^ yVa.ln^-)" / dz 1 dz 2 ...dz n 

x A n (x, z\ z 2 , ...z n , </)Tr{t/£° [x, l}-Ul Cn [1, x]+} 

x Tr{f/f° [1, 2]_f/ 2 tCo [2, 1]+} . . . Tr{U$> [n, y]-U^ [y, n} + } 

+ B n {x,Z X ,Z 2 ,...Z n ,y) 

x Tr{C/i° [x, l]-Ul C ° [1, 2]+C/ 2 Co [2, 3]_C/ 3 tCo [3, l] + E/f° [1, 2]_(7 2 tCo [2, x}+} 
x Tr{C/ 3 Co [3,4]_C/ 4 tCo [4,3] + }...Tr{C/^Ky]_^[ 2/ ,n] + } + ... 



+ N^C n (x,z\z 2 



\y;)Tr{U^[x,y}.U^[y,x} + } 
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(153) 



where U$ ] = C/ (t) (^), = [xi,Xj] and 

A n (x, z , z 2 , . . . , z n , y), B n (x, Z , z 2 , . . . , z™, y), . . ., C n (x, z , z 2 , . . . , z" , y) a re 
the meromorphic functions that can be obtained by using the Eqs. (|300| , |30l[ ) n 
times which give us a sort of Feynman rules for calculation of these coefficient 



functions. If we now evol ve ou r operators from £ ~ to given by Eq. (143) 
we shall obtain a series fll53|) of matrix elements of the operators (U) n (W) n 
normalized at Co- These matrix elements correspond to small energy ~ m 2 and 
they can be calculated either perturbatively (in the case the "virtual photon" 
matrix element ) or using some model calculations such as QCD sum rules in 
the case of nucleon matrix element corresponding to small- x 7*p DIS . It should 
be mentioned that in the case of virtual photon scattering considered above 
we can calculate the matrix elements of operators l/W ...l/W perturbatively. 
Because U = 1 + ig jA^dx^ + ..., in the leading order in a s we can replace by 1 
all but two U(W)'s, so we return to the BFKL picture describing the evolution 



of the two operators l/W . The non-linear equation (137) enters the game in 
the situation like small-x DIS from a nucleon or nucleus when the matrix 
elements of the operators UW ...l/W are non-perturbative, consequently there 
is no reason to expect that extra U and U> will lead to extra smallness. In 
this case, at the low "normalization point" Co one must take into account the 
whole series of the operators in the right-hand side of Eq. ( |l53| ), indicating 
the need for all the coefficients a ni b n ...c n . Recently, these coefficients were 
calculated by Y. Kovchegov E3 for the case of DIS from the large nuclei in 
the McLer ran- Venugopalan model, and the results indicate that th£ non-linear 



equation (137) leads to unitarization of the pomeron in this caseEZI 



The zoo of different Wilson-line operators ( |153| ) may be reduceded by using 



the dipole picture Ol3 Technically, it arises when in each order in a s ln(^) we 

keep only the term Tr{U^°i7| Co }Tr{i7f o U'] Co }...Tr{Z7^L^ C[ '}-subtractions0 in 
right-hand side of Eq. ( |153[ ); for example, in Eq. ( |152| ) we keep the two first 
terms and disregard the third one. In other words, we take into account only 
those diagrams in Fig. [59| which connect the Wilson lines belonging to the 
same Ti{UkUl +1 }. (This corresponds to the virtual photon wave function 
in the large- N c approximation). The diagrams of the corresponding effective 



theory are obtained by multiple iteration of Eq. (137) and give a picture where 
each "dipole" Tr{[/fcl^ +1 } can create two dipoles according to Eq. (137). The 
motivation of this approximation is given in Refs. 24, 25, and the discussion 
of unitarization of the BFKL pomeron in the dipole picture is presented in 
Ref. 41. 



m By "subtractions" we mean this operator with some of the Tr{[/ fc l/t +i } substituted by 
N c . 
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3.7 Operator expansion for diffractive high-energy scattering 



The nonlinear term in the equation (137) describes the triple vertex of hard 



pomerons in QCD. In order to see that, it is convenient to consider some 
process which is dominated by the three-pomeron vertex — the best example 
is the diffractive dissociation of the virtual photon. 

The relevant operator expansion for diffractive scattering is .obtained by 
direct generalization of our approach to the diffractive processesEa The total 
cross section for diffractive scattering has the form: 



riiff 



d ^ qx J tJ^EWmW + x W + *l>(o)b>, (154) 



p and p' are the nucleon momenta and J^x means the summation over all 
the intermediate states. We can formally write down this cross section as a 
"diffractive matrix element" (cf. Ref. 43): 



^ = W^= I dxe^(p\T{j-(x)j+(0)}\p), (155) 



where [] 



def 



(p\T{3-3+ef dz(C+ ^~ C -^}\p) (156) 

x (p' + X\T{M0yf dzC(z) }\p). 

The superscript "-" marks the fields to the left of the cut and + to the right. 
The definition of the T-product of the fields with ± labels is as follows: the 
+ fields are time-ordered, the — fields stand in inverse time order (since they 
correspond to the complex conjugate amplitude), and — fields stand always 
to the left of the + ones. Therefore, the diagram technique with the double 
set of fields is the following: contraction of two + fields is the usual Feynman 
propagator p2 ^_ ie (for the quark field), contraction of two — fields is the com- 
plex conjugated propagator $_ i€ , and the contraction of the — field with the 
+ one is the "cut propagator" 2ir6(p 2 )9(po) This diagram technique for 
calculating T-products of double set of fields exactly reproduces the Cutkosky 



"The differen ce b etween Eq. (154}and the last line in Eq. (1159) is that j's sue Heisenberg 



operators in (154) while in Eq. (156) the operators stand in the interaction representation 
"We will use the — h perturbative propagator only for hard momenta, hence the additional 
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rules for calculation of cross sections. The light-cone expansion of the diffrac- 
tive matrix element (155) gives operator definition of the diffractive parton 
distributions^ 

Let us discuss the high-energy operator expansion for the diffractive am- 
plitude W^jf. Similarly to the case of usual amplitude (110), we get in the 
lowest order in a s 

flavors 

x {NlTriW^^i^W^'^/'i-kx^lN), (158) 

where W(k±) is a Fourier transform of 

W(x ± ) = V\x ± )U(x ± ), W*(xx) = Ut(xx)V(x x ). (159) 

Here U(x±) denotes the Wilson- line operator constructed from + fields and 
V(x±) denotes the same operator constructed from — fields: 

U^(xj_) = [oopx +xj_, — oopi + , V(xj_) = [ooe + a;_L, — ooe + sjj - . (160) 

After integration over fast quarks, the slope of the Wilson lines is £ = xb = 
see the discussion in Sec. 3.3. 

The evolution equation (with respect to the slope of the supporting line) 
turns out to have the same form as Eq. (137) for non-diffractive amplitudes: 



/- d n a ° N c f j {$-y)\ ,,,,, 

C T( w(x^) =-^rj dzx {S _^_ yli (i6i) 

x {W{x ± , zi _)+W(x x ,z x ) - W{x ± ,y x ) + W{x u z 1 _)W(z u y x )}, 
where 

W(x x ,y x ) d = ^Tr{W(x ± )W^(yx)} - 1, (162) 



emitted nucleon with momentum p' (constructed from soft quarks) can be factorized 

W(x)\p' + X)(p' + X|^(0)|0> (157) 

X 

~ W(*)|x}<x|^(o)|q) ® \ P ')(p'\ = J (U^t ^<5(p 2 )e( P „) ® \ P ')( P '\. 
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(cf. Eq. (138). Similarly, the linear evolution is: 



(N\W Cl (x ± ,0)\N) 



dv 
2^ 



(163) 



x dz ± (zl)-?-™(N\W^(z ± ,0)\N}, 



where oj(v) — 2^-N c x(v), see Eq. (|56|). Let us now describe the diffractive 
amplitude in LLA and in leading order in N c . In th is approximation we must 
take into account the non-linearity in the Eq. ( 161 ) only once, the rest of the 
evolution is linear. The result is (roughly speaking) the three two-gluon BFKL 
ladders which couple in a certain point, see Fig. EOl For the case of diffractive 





P p 

Figure 20: Amplitude of diffractive scattering in the LLA-7V C approximation. 

DIS, this evolution has the form (cf. Ref. 45): 

(N\ [ d yi _W c = XB (x x +y±,y±)\N) (164) 



^ f dvdv 1 dx 1 dv 2 dx 2 {x 2 L )^ +w ((x! - ^ly^^ 1 ^ 2 ^ 
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f d 2 p' ± r dM 2 t s \<"W fM 2 



#2 , „ /,/2\ a ("l)+«h) 



TV 

x (p\U^(x 1 ^ 1 )\p'){ P '\U^(x 2 ,v 2 )\p), 
where M 2 is the invariant mass of the produced particles, and 

r(| - i{v + vi- V2))T{\ - i(v + u 2 )) 



6(i/; ui, v 2 ) 



r(i + i(y + v x - v 2 ))Y{\ + i(y + u % )) 



x T ™kV v \ ^\ +iv '\- ivi '\- iv ^ (165) 

I + iv) 2 I 2 

is a certain numerical function of three conformal weights (the explicit form 
was found in Ref. 46 ) which has a maximum 6(0, 0, 0) = 2tt 7 4F3(^) 6 F 5 (^) ~ 
7766.679 . The value of M 2 determines the rapidity gap: from rj = In to 

X] = In we have a production of particles described by the cut part of the 

ladder in Fig. ^ which brings in the factor [s j M 2 ^ 1 ^ while from rj = In 
to rj = \tixb we have a rapidity gap so there are two independent BFKL 
ladders which bring in the factors (M 2 /Q 2 ) uj( " i) and (M 2 /Q 2 ) u(v2) . Since the 
intercept of the BFKL pomeron ojq > 0, this cross section increases with the 
growth of the rapidity gap. 

The coupling of BFKL ladder with non-zero momentum transfer to a nu- 
cleoli is described by the matrix element (p'\U(x, v)\p). As we discussed in the 
previous section, at large momentum transfer it can be approximated by the 
electromagnetic form factor of the nucleon, 

(g - y 1 )! 

2ms / 
If one interpolates the form factors by the dipole formulas, the diffractive 



w(*,„) = l^dy'[ {sl _ m ^_ s? J prr^jr (166) 



amplitude in the LLA-7V C approximation (164) can be calculated numerically. 



The non-linar equation (161) can be applied to the diffractive DIS from 
the nuclei. In this case there is an additional large parameter, the atomic 
number A, and therefore one should take into account the multitude of the 
non-linear vertices rather than one vertex as in Fig. These "fan" diagrams 
were summed up in Ref. 47 resulting in a cross section which has a maximum 
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at a certain rapidity gap (unlike the ~LLA-N C model for the nucleon where the 
cross section increases with the rapidity). 



4 Factorization and effective action for high-energy scattering 

4-1 Factorization formula for high-energy scattering 

Unlike usual factorization, the coefficient functions and matrix elements enter 
the expansion ( |80| ) on equal footing. We could have integrated first over slow 
fields (having the rapidities close to that of ps) and the expansion would have 
the form: 

A M =X! / rf 2 ^i---rf 2 ^^ 1 --""(^i,--^n)(p J 4|Tr{i7 Jl ( a ; 1 )...;7 4?i ( a; „)}|^). 

(167) 

In this case, the coefficient functions D are the results of integration over 
slow fields ant the matrix elements of the U operators contain only the large 
rapidities r/ > 770. The symmetry between Eqs. (1) and (2) calls for a factor- 
ization formula which would have this symmetry between slow and fast fields 
in explicit form. 

I will demonstrate thatxine can combine the operator expansions (^0|) and 



( 167 ) in the following wayEa 



A(s,t) = J d 2 xi...d 2 x n (168) 

x ( P A|^ ai<1 (a:i)...^ a » i -(a;„)|y A )< P fl|^Hx 1 )...tC'(a:«)|l4), 

where Uf = Ti(X a Ui) (A a are the Gell-Mann matrices). It is possible to rewrite 
this factorization formula in a more visual form if we agree that operators U act 
only on states B and B' and introduce the notation Vi for the same operator 
as Ui only acting on the A and A' states: 

A(s,t) = (p A \(p B \exp(iJd 2 xV ai (x)U?(x)) \p' A )\ P ' B ). (169) 

The supporting lines of both U and V operators are collinear to the vector 
n corresponding to the "rapidity divide" rj . The explicit form of this vector 



is n = upi + crp2, where a = — and ]na/a = r\. In a sense, formula (169) 



amounts to writing the coefficient functions in Eq. ( |80| ) (or Eq. (167)) as matrix 
elements of Wilson-line operators. Eq. ( |169| ) illustrated in Fig. |2l| is our main 
tool for factorizing in rapidity space. 

In order to understand how this expansion can be generated by the factor- 



ization formula of Eq. (169) type we have to rederive the operator expansion in 
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Figure 21: Structure of the factorization formula. The vector n gives the direction of the 
"rapidity divide" between fast and slow fields. 



axial gauge A, — with an additional condition A*| = _ co = (the existence 
of such a gauge was illustrated in Ref. 49 by an explicit construction). It is 
important to note that with with power accuracy (up to corrections ~ a) our 
gauge condition may be replaced by n^A^ — 0. In this gauge the coefficient 
functions are given by Feynman diagrams in the external field 

Bi{x) = Ui(x±)Q(x,), B. = B,=0, (170) 

which is a gauge rotation of our shock wave (it is easy to see that the only 
nonzero component of the field strength tensor F,i{x) = Uj{x±_) 8{x* ) corre- 



sponds to shock wave). The Green functions in external field (170) can be 
obtained from a generating functional with a source responsible for this exter- 
nal field. Normally, the source for given external field A^ is just J„ = D^F^, 
so in our case the only non- vanishing contribution is J*(i?) = D l F^. However, 
we have a problem because the field which we try to create by this source does 
not decrease at infinity. To illustrate the problem, suppose that we use another 
light-like gauge A* — for a calculation of the propagators in the external field 



(170). In this case, the only would-be nonzero contribution to the source term 
in the functional integral D l F,A* vanishes, and it looks like we do not need 
a source at all to generate the field (This is of course wrong since B^ is 
not a classical solution). What it really means is that the source in this case 
lies entirely at the infinity. Indeed, when we are trying to make an external 
field A in the functional integral by the source J M we need to make a shift 
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An — > + .Au in the functional integral 



/ VA exp - i/d 4 a;J°(a;)^(x)} 



(171) 



after which the linear term D p, F iLV A v cancels out with our source term J^A^ 
and the quadratic terms lead to the Green functions in the extern al fie ld A. 
(Note that the classical action S(A) for our external field A — B ( ]170| ) van- 
ishes). However, in order to reduce the linear term J d xF^D^Av in the 
functional integral to the form J d^xF) lL F jLV A v (x) we need to perform an in- 
tegration by parts, and if the external field does not decrease there will be 
additional surface terms at infinity. In our case we are trying to make the 
external field A = B, consequently the linear term which need to be canceled 
by the source is 



/ dx,dx*d 2 x±Fi.D*A % — J dx*d 2 x_LF l ,A t |^" = !° c 



(172) 



This contribution comes entirely from the boundaries of integration. If we 
recall that in our case F,i(x) = Ui(x±)S(x*) we can finally rewrite the linear 
term as 



/ d 2 x±U i (x±){A i (-oop 2 + x±) - A i (oop 2 + xx)}- 



(173) 



The source term which we must add to the exponent in t he fu nctional integral 
to cancel the linear term after the shift is given by Eq. (173) with the minus 
sign. Thus, Feynman diagrams in the external field (170) in the light-like gauge 
A* = are generated by the functional integral 

J VAexp {iS(A) + % J d 2 x ± U ai {x ± )[AU°°P2 + x ± ) - A ai {-oo P2 + xl)]}. 

(174) 

In an arbitrary gauge the source term in the exponent in Eq. (174) can be 
rewritten in the form 

P /"OO 

2i / d 2 x±Tr{U z (x±) / dv[-oop2,vp 2 }x ± F^{vp2 + x±)[vp 2 ,~oop2]x ± }- 



(175) 

Therefore, we have found the generating functional for our Feynman diagrams 
in the external field (170). 

It is instructive to see how the source (175) creates the field ( |17C| ) in 
perturbation theory. To this end, we must calculate the field 



A^x) 



VAA^x) exp \iS(A) + 2iJ d 2 x ± Tr{U l {x±) 

dv[-ocp 2 ,vp2] x _ L F*i(vp2 + x±)[vp2, -oop 2 W}} (176) 



59 



by expansion of both S(A) and gauge links in the source term (175) in powers 
of g (see Fig. |22|). In the first order one gets 





(c) 



Figure 22: Perturbative diagrams for the classical field (170) 



A^(x) =i dv dz ± U ia (z ± )(A^x)F.^(v P2 + z±)), (177) 

J — oo J 

where (O) = jVAe' lSo O. Now we must choose a proper gauge for our cal- 



culation. We are trying to create a field (17C) perturbatively and therefore 



the gauge for our perturbative calculation must be compatible with the form 



(170), otherwise, we will end up with the gauge rotation of the field B{x). (For 
example, in Feynman gauge we will get the field A^ of the form of the shock 
wave Ai = A* = 0, A, ~ 5(x*)). It is convenient to choose the temporal 
gauge Ao — ^| with the boundary condition •A\ t __ OB = where 

A^x) = f dt'F 0ll (t',x). (178) 



p The gauge At = which we used above is too singu lar fo r the perturbative calculation. 
In this gauge one must first regulate the external field ( |l70[ ) by the replacement Ui9(x t ,) — > 
Ui9(x t )e~ ex ' and let e — » only in the final results. 
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In this gauge we obtain 

7(0)/ \ _ f d P ( Py(Pl +P2)u + (P V) . ^VpPu \ 

^ [X) ~ J (2nr s ( a + /3 + ie) + s(a + /? + itf ) 



1 



. dz ± e- iax '- i P x * +i ^-^p 2 J(a^)d i U ja (z 1 _) (179) 
aps — pj_+ te J 2 

where <5(a|) comes from the J dv e m a % . (Note that the form of the singularity 
( Po +ie) which follows from Eq. ( |178|) differs from the conventional prescription 
V.p.-^). Recalling that in terms of Sudakov variables dp = ^dad(3dp± one 
easily gets ^1°^ = *4i°^ = and 

Af\x) = 0(x.) J J^^r J dz ± e^-^d i d j U^{z±), (180) 



or more formally, 

= -0(t~)- 
dl 

= U i (xA.)9(x,,)-9{x*)\0 l x gi i + did j )U i (xi_\ (181) 
d\ 

(in our notations d\ = —did 1 ). Now, since Ui(x) is a pure gauge field (with 
respect to transverse coordinates) we have diUj — djUi — i[Ui, Uj], so 

A<f\x) = Ui{x±)e(x*) - 6(x*)ig^[Ui, Uj])(x±). (182) 



Consequently, we have reproduced the field (|l70j) up to the correction of of g. 
We will demonstrate now that this 0{g) correction is canceled by the ne xt-to - 
leading term in the expansion of the exponent of the source term in Eq. ( |l 76| ) . 
In the next-to-leading order one gets (see Fig. p2|b) 



A^(x) = g jdy Jdz ± dz' ± W a (z ± )U kb (z' ± ) (183) 
x (a, {x)2Tt{d a A P (y) [A a (y),Ap(y)}} 
x { dvF^(v P2 + z±) j dv'F» k (vp 2 + z' ± )). 
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It is easy to see that „4i^ = A^ = and 



Since 4°^ is given by Eq. ( |l82|) , this reduces to 
A { l\x) = ~ g e(x*)jdy ± ^ ^ i&mMMy)]) + 0(g 2 ). (185) 



The right-hand side of this expressions cancels the second term in Eq. (182) 
and we obtain 

Ai(x) = Ui(x±)0(x.) + 0(g 2 ). (186) 

Likewise, one can check that the contributions ~ g 2 coming the diagrams in 
Fig. p2]c cancel the g 2 term in the Eq. (186). Taking into account arbitrary 



number of the tree-gluon vertex iterations, one gets the expression Ui(x±_)9(x*) 
without any corrections. 

We have found the generating functional for the diagrams in the external 



field (170) which give the coefficient functions in front of our Wilson-line oper- 
ators fTj. Note that formally we obtained the source term with the gauge link 
ordered along the light-like line, a potentially dangerous situation. Indeed, it 
it is easy to see that already the first loop diagram shown in Fig. ^3] is diver- 
gent. The reason is that the longitudinal integrals over a p are unrestricted 
from below (if the Wilson line is light-like). However, this is not what we want 
for the coefficient functions because they should include only the integration 
over the region a p > a (the region a p < a belongs to matrix elements, see 
the discussion in Sec. 3). Therefore, we must impose somehow this condition 
a p > a in our Feynman diagrams created by the source (|l75| ). Fortunately 
we already faced similar problem — how to impose a condition a p < a on the 
matrix elements of operators U (see Fig. |l5|) — and we solved that problem 
by changing the slope of the supporting line. We demonstrated that in order 
to cut the integration over large a > a from matrix elements of Wilson-line 
operators Ui we need to change the slope of these Wilson-line operators to 
n = api + ap2- Similarly, if we want to cut the integration over small a p < a 



from the coefficient functions we need to order the gauge factors in Eq. (175) 
along (the same) vector n = api + <7p 2 -n 



*?Note that the diagram in Fig. |23| is the diagram in Fig. |15| b turned upside down. In the 
Fig. hab diagram we have a restriction a < a. It is easy to see that this implies a restriction 
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] liglire 23: A typical loop diagram in the external field created by the Wilson-line source 
(175). 



Therefore, the final form of the generating functional for the Feynman 



diagrams (with a > a cutoff) in the external field (170) is 

J VAD^cxp ^iS{A, *) + i J d 2 x x U ai (x ± )V^(x x )\ , (187) 

where 

Vi(x±) = / dv[— oon, vn] x n^F^i{vn + x±)[vn, — oon] x , (188) 



and V°" = Tr(A a Vi) as usual. For completeness, we have added integration 
over quark fields so S(A, is the full QCD action. 

Now we can assemble the different parts of the factorization formula ( |169| ) . 
We have written down the generating functional integral for the diagrams 
with a > a in the external fields with a < <r; what remains now is to write 
down the integral over these "external" fields. Since this integral is completely 



independent of (187) we will use a different notation B and x f° r t ne a < a 
fields. We have 

VAV^e iS ^j{p A )j{p' A )j{p B )j{jp' B ) (189) 



j3 > <T if one chooses to write down the rapidity integrals in terms of /3's rather than a's. 
Turning the diagram upside down amounts to interchange of and pg, leading to (i) 
replacement of the slope of the Wilson line by api + ap2 and (ii) replacement a <-> f3 in the 
integrals. Thus, the restriction /3 > a imposed by the line collinear to ap\ + ay ? in diagram 
in Fig. 15b means the restriction a > a by the line || ap\ + ap2 in the Fig. E3 diagram. 
After renaming cr by <r we obtain the desired result. 
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VAVi>V^ s{A ^ j ( PA ) j (p' A ) / VBVxDx 
3(PB)3(p B y SiBa) exp {i ( d 2 x ± U ai (x x )V?(x ± )}. 



The operator Ui in an arbitrary gauge is given by the same formula ( |l88| ) as op- 
erator Vi with the only difference that the gauge links a nd F ti are constructed 



from the fields B^. This is our factorization formula (169) in the functional 
integral representation. 

The functional integrals over A fields give logarithms of the type g 2 In 1/a 
while the integrals over slow B fields give powers of g 2 ln(as/m 2 ). With loga- 
rithmic accuracy, they add up to g 2 In s/m 2 . However, there will be additional 
terms ~ g 2 due to mismatch coming from the region of integration near the 
dividing point a ~ a, where the details of the cutoff in the matrix elements of 
the operators U and V become important. Therefore, one should expect the 
corrections of order of g 2 to the effective action J dx±U i Vi of the type 

exp (i / d 2 x±Ui(x±_)Vi(x±) + i / dx±dy±dz± (190) 



x 



Ui(x±)Ui{y±)Vi(z±)Vi(t±)K(x± -t±,y±-t±,g±-t±)j 



where K is a calculable kernel. In general, the fact that the fast quark moves 
along the straight line has nothing to do with perturbation theory (cf. Ref. 50) , 
therefore it is natural to expect the non-perturbative generalization of the 
factorization formula constructed from the same Wilson-line operators Ui and 
Vi. 

4-2 Effective action for given interval of rapidities 

The factorization formula gives us a starting point for a new approach to the 
analysis of the high-energy effective action. Consider another rapidity n' in 
the region between 770 and t\b = In m 2 / s. If we use the factorization formula 



(|189|) once more, this time dividing between the rapidities greater and smaller 



than n' , we get the expression for the amplitude (|6j) in the form (see Fig. |2J): r 

iA(s,t) = JvAe* s ^j(pA)](p' A )](PB)j(p'B) (191) 
= [vAe iS Wj(p A )j(p' A ) JvBe^^pB^ip's) 



""Strictly speaking, the l.h.s. of Eq. (191) contains an extra 16-7r 4 <5(p j 4 + p' A — ps — p' B ) in 
comparison to the amplitude (ft). 
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x J VCe tS< - c ^e i J d2x ± vai ( x ±) Y i a ( x ±) +i J d 2 x ± w ai (x ± )u?(x ± ) 
(For brevity, we do not display the quark fields.) In this formula the operators 




Figure 24: The effective action for the interval of rapidities rjo > r\ > rj' . The two vectors n 
and n' correspond to "rapidity divides" rjo and rj' bordering our chosen region of rapidities. 



Vi (made from A fields) are given by Eq. (188), the operators Yi are also given 
by Eq. fll8S| ) but constructed from the C fields instead, and the operators Wi 
(made from C fields) and Ui (made from B fields) are aligned along the direction 
n' = a'pi + a'p2 corresponding to the rapidity 77' (as usual, Inc'/cr' = 77' where 
a' = m 2 /sa'), 



Vi{A) x± = j dv[-aon, vn] x n li F li i(vn + x±)[vn, -oan] x , 
J —00 

/oo 
dv[-oon,vn] x n fJ 'F fli (vn + x±)[vn, —oon] x , 
-00 

/oo 
dv[-oon', vn'] x n T^vn' + x±)[vn', -oon'] x , 
-00 

/•OO 

Ui(B) x± = / dv[-oon',vn'] x n ^F^vn' + x±)[vn', -oon%. 



In conclusion, we have factorized the functional integral over "old" B fields 
into the product of two integrals over C and "new" B fields. 



65 



Now, let us integrate over the C fields and write down the result in terms 
of an effective action. Formally, one obtains: 



iA(s,t) = JVAe^j(p A )j(p' A ) jVBe^]{p B )j{p' B y s ^ u ^\ (192) 
where the effective action for the rapidity interval between 77 and rj' is defined 



as 



(Ui = Wi-diU and Vi = V'~diV as usually). This formula gives a rigorous 
definition for the effective action for a given interval in rapidity. 

Next step would be to perform explicitly the integrations over the longi- 



tudinal momenta in the right-hand side of Eq. (193) and obtain the answer 
for the integration over our rapidity region (from t/q to t]q) in terms of two- 
dimensional theory in the transverse coordinate space,^] hopefully resulting in 
the unitarization of the BFKL pomeron. At present, the known how to do this. 
One can obtain, however, a first few terms in the expansion of effective action 
in powers of Vi and Ui. The easiest way t o do this is to expand gauge factors 
Yi and Wi in right-hand side of Eq. (|l93| ) in powers of C fields and calculate 



the relevant perturbative diagrams (see Fig. 25). The first few terms in the 



(a) 



(b) 



(c) 



(d) 



Figure 25: Lowest order terms in the perturbative expansion of the effective action. 

effective action at the one-log level | have the formij^l 

s Historically, the idea how to reduce QCD at high energies to the two-d imen sional effective 



theory was first suggested in Ref. 5f where the leading term in Eq. (194) was obtained. 
However, careful analysis of the assumptions made in this paper shows that the authors 
considered the fixed-angle limit of the theory (s, t — > 00) rather than the Regge limit (where 
— > 00 but t is fixed), ft turns out that the first term in Eq. (194) is the same for both limits, 
but the subsequent terms differ. 

*This "one-log" level corresponds to one-loop level for usual Feynman diagrams. Supcrfi- 
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(a) (b) (c) (d) 

Figure 26: Counting of loops for Feynman diagrams (a),(c) and the corresponding Wilson- 
line operators (b),(d). 



S cS = J d z xV ai (x)U?{x) (194) 
In - f [N c [ d 2 xd 2 yV^(x) \n 2 (x ~ yfU^y) 



647r 3 a 

fabcfi 



t j d i xd?yd*x>Wd 2 zVZ i (x)V%(y)Ul k tf)UZ l tf) 
(x - z) 2 ^ (y- zf (d \ 2 ^ (x> - z) 2 ^ (y> - z) 2 \ 



(x — x') 2 (y — y') 2 \dzi) (x — x') 2 (y ~ y') 2 I 

where we we use the notation V?Ax) = ■^-V i a (x) etc. The first term (see 
Fig. |25|a) looks like the corresponding term in the factorization formula ( |189| ) , 
only the directions of the supporting lines are now strongly different |] The 
second term shown in Fig. E5c is the first-order expression for the reggeization 



cially, the diagram in Fig. |25p looks like tree diagram in comparison to diagram in Fig. |25p 
which has one loop. However, both of the diagrams in Fig. and d contain integration over 
longitudinal momenta (and thus the factor In ) so in the longituduinal space the diagram 
in Fig. pj| d is also a loop diagram. This happens because for diagrams with Wilson-line 
operators the counting of number of loops literally corresponds to the counting of the num- 
ber of loop integrals only for the transverse momenta. For the longitudinal variables, the 
diagrams which look like trees may contain logarithmical loop integrations. This property 
is illustrated in Fig. the Wilson-line diagram shown in Fig. pjj b has two loops and the 
diagram shown in Fig. pi| d is a tree but both of them originated^ from Feynman diagrams 
shown in Fig. 26b and c with equal number of loops. To avoid confusion, we will use the 
term "one-log level" instead of "one-loop level." 

"Strictly speaking, the contribution coming from the diagram shown in Fig. baa has the form 
J d 2 xV at (x) d 'J^ j U a: > (x) which differs from the first term in the right-hand side of Eq. (194) 
by J d 2 xV ai {x)-^(d 2 g l j - d,dj)U a Hx). Yet, it may be demonstrated that this discrepancy 
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of the gluon (|74j) and the third term (see Fig. |25]d) is the gluon emission term 
in the BFKL kernel ( p6| ) in the impact parameter representation. 

Let us discuss subsequent terms in the perturbative expansion (194). There 
can be two types of the logarithmical contributions. First is the "true" loop 
contribution coming from the diagrams of the Fig. p7| a type. This diagram is 
an iteration of the Lipatov's Hamiltonian. In addition, in the same (In -p-) 2 or- 
der there is another contribution coming from the diagram shown in Fig. |27|b. 
In perturbation theory, these two contributions are of the same order of mag- 



(a) (b) 
Figure 27: Typical perturbative diagrams in the next (in order. 



nitude. 

The situation is different for the case of scattering of two heavy nuclei. 
Assuming that the effective coupling constant is still small due to the high 
densityp we see that jCl, yet the sources are strong (~ |) so gUt ~ gYi ~ 1. 

In this case, the diagram in Fig. |27| a has the order g^U^V 2 (in p-) 2 ~ (in p) 2 



while the "tree" Fig. 25 b diagram is 

~ 9 4 U?V? (m^) 2 ~ l(lr^) 5 (195) 



In this approximation, first we shall sum up the tree diagrams. As usual, the 
best way do this is to use the semiclassical method which will be discussed in 
Sec. 5. In the next paragraph we will consider the intermediate situation with 
one weak source and one strong source. 



(which is actually ~ O(g) for a a pure gauge field Ui) is canceled by the contribution from 
the diagram with the three-gluon vertex shown in Fig. pj| b just as in the case of perturbative 
calculation of Ai discussed in Sec. 3. 
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4-3 Effective action for one weak and one strong source 



Consider again the DIS from a nucleon or nucleus where the high-energy be- 
havior is governed by the non-linear evolution equation (137). In this section 
we will translate the evolution results (153) into the effective action language 
(see also Refs. 52, 53). In the case DIS one of the sources (corresponding to 
quark-antiquark pair) is weak while the other (describing the nucleon or nuclei) 
is strong. 

For example, if the source Vi is weak (and hence gVi is a valid small 
parameter) but the source Ui is not weak (so that gVi ~ 1 is not a small 
parameter), one must take into account the diagrams shown in Fig. p8| a and 
b. The multiple rescatterings in Fig. p8|a,b describe the motion of the gluon 



(a) (b) 

Figure 28: Perturbative diagrams for the effective action in the case of one weak source and 
one strong one. 

emitted by the weak source V* in the strong external field Ai = Ui9{x*) created 
by the source Ui. The result of the calculation of the diagram in Fig. |2§| a 
presented in a form of the evolution of the Wilson-line operators Ui can be 



easily obtained using the evolution equations (301) 

U?(x±) - U?(x ± ) (196) 

- & ln 7? J dy^ w 1 ^(f abc mu y ) bc + N c un^)) + ---, 

where dots stand for the terms with higher powers of g 2 \n-p. This evolution 
equation means that if we integrate over the rapidities t/q > r\ > t/q in the ma- 



trix elements of the operator Yi we will get the expression (196) constructed 



from the operators Ui with rapidities up to r)' times factors proportional to 
g 2 (ijo — n' ) = g 2 In p.. Therefore, the corresponding contribution to the effec- 
tive action at the one-log level takes the form 

[ dx ± V?(x ± )U ai (x x ) -» / dx±V i a (x±)U tti (x ± ) (197) 
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£-7 In — I dx±dy± 

S7T J & 1 



{y i {x 1 _)Utd i U y ) aa - N c V a \x ± )U?(x ± )) 



where the first term is the lowest-order effective action (= the first term in 
Eq. (194)) and the second term contains new information. To check the second 
term, we may expand it in powers of the source Ui, then it is easy to see that 
the first nontrivial term in this expansion coincides with the gluon-reggeization 
term in Eq. (194) . 

Apart from the (197) term, there is another contribution to the one-loop 
evolution equations coming from the diagrams in Fig. p8|b. It can be easily 
obtained using formulas (30C) from the Appendix, 



U?(x±)Ufoj 



9 , 

47T 3 & 



dz± 



(x- z,y- z) 
(x-z) 2 x (y- z) 



(198) 

y \ab 



*r{ulu y + i-ulu z -ulu y )] V,-) , 



where 



VfO(x ± ) 
0{y±_) v! 



jO(x ± )-iUi( Xl _)0{x±), 



d_ 

dy" 



■0(y ± )-iO(y ± )Ui{y ± ), 



(199) 



are the "covariant derivatives" (in the adjoint representation). The correspond- 
ing term in effective action is 



^ln- / dx ± dy ± (Vf^°) (x ± ) / d Zl _ 



(x- z,y- z)j_ 
{x-z)\{y- zf 



x (UlU y + l-UtU z -UtU y ) ab (^)(y ± ) 



(200) 



The final form of the one-log effective action for this case is the sum of the 
expressions (197) and (|200|) , 



S${Vi, Uj) = J d 2 xV ai (x)U?(x) + ln£ J dx ± dy ± j^- 

{i(V\x 1 _)Uld % U y ) aa - N c V ai (x ± )U?(xx)) 

(x-z)± ■ (y- z)± 



— r In — 

8tt 3 a' 



dx ± dy ± V*V a \x ± ) / dz ± 



(x-z)±(y- zf_ 



x (Uiu y + 1 - V$U, UtU y ) a \nv b i{y x ), 



(201) 
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where Vi is a weak source and Ui is a strong one. It is clear that if the source 
Vi is strong and Ui is weak diagrams the effective action S^'(Vi, Uj) will have 
the similar form with the replacement V <-> U coming from the diagram shown 
in Fig. M. 



Figure 29: Effective action for the strong source V and the weak source U. 

As we mentioned above, in the case of two strong sources the (in -p) 



terms start from the diagram shown in Fig. |27| b (see Eq. ( |195| )), hence Fig. 27 
and Fig. |2^ complete the list of diagrams which contribute to the effective 
action at the one-log level. Higher-order diagrams start from higher powers 
of In -2? . The analog of LL A here is a cluster expansion with the parameter 
(U — 1)(V — 1) In p shown in Fig. |3(j. Of course, the diagrams of Fig. ^ give 



Figure 30: Cluster expansion of the effective action. 

the terms ~ In p too, but in the leading order the kernel of the corresponding 
evolution equation is determined by Fig. ^ and Fig. ^8|. Thus, the one-log 
answer for two strong sources can be guessed by comparison of the answers 
for S eS (Vi,Uj) with V ~ 1, U l - | and with Ui ~ 1, Vi ~ |. Instead of 
doing that, we will obtain the one- log result for two strong sources using the 



semiclassical method and check that it agrees with (201). 



It means that the one-log answer in the general case can be guessed by 
comparison of the answers for S e s(Vi, Uj) with Vi ~ 1, Ui ~ ^ and with 
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Ui ~ 1, VI ~ i Instead of doing that, we will obtain the one-log result for two 
s tron g sources using the semiclassical method and check that it agrees with 



5 High-energy effective action in sQCD 

5. 1 Effective action and collision of two shock waves 



The functional integral (193) which defines the effective action is the usual 
QCD functional integral_with two sources corresponding to the two colliding 
shock waves, see Fig. UWj Instead of calculation of perturbative diagrams we 



E»m 




Shockwaves 



Figure 31: Scattering of two shock waves. 



can use the semiclassical approach which is relevant when the coupling constant 
is relatively small but the characteristic fields are large - in other words, when 
g 2 <C 1 but gVi ~ gUi ~ 1. As was discussed in Ref. 4, this situation is 
realized in the heavy-ion collisions where the coupling constant is defined by 
the parton saturation scale Q s , which is estimated to be ~ 1 GeV at RHIC and 
^2 — 3 GeV at LHCH0 Even if we consider the 7*7* scattering, the number 
of gluons in the middle of the rapidity region may become very large leading 
to the saturation at high energies so in the middle of the rapidity region we 
will se the scattering of two strong shock waves. 



If both sources are strong, one can calculate the functional integral (193) 
by expansion around the new stationary point corresponding to the classical 
wave created by the collision of the shock waves. With leading log accuracy, 
we can replace the vector n by p\ and the vector n' by P2 ■ Then the functional 



integral (193) takes the form 



JS et t{V,U; 



(202) 
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where now 

/OO pOG 

dvF 9i (vp 1 + x±), W?= dvF«{vpa+x±). (203) 
- oo J —oo 

Hereafter we use the notations 

0(x) = [— oapi + x, x]0(x)[x, — oapi + x], 

6{x) = [-oop 2 + x,x]0(x)[x, ~oop 2 + x]. (204) 

Note that we changed the name for the gluon fields in the integrand from C 
back to A. 

As usual, the classical equation for the saddle point A in the functional 
integral ( p02] ) is 

5 



SA 



Sqcd + fd 2 x±V ai (x±)Y i a (x±) + J d 2 x±_W ai U?(x±_ 



= 0. (205) 

A=A 



To write them down explicitly we need the first variational derivatives of the 
source terms with respect to gauge field. We have: 



SYi = SAi(oopi + x±) — 5Ai(—oopi + x±) — / duV i5Ai(upi + x±), 

J —oo 

/oo 
duViSAi(up2 + x±), (206) 
-oo 

where 

V t O(x) = diO{x)-i{Y i {xi_)+A i {-oopi + x 1 _),0{x)], 

fiO(x) ee d i O(x)-i[W i (xx)+A i (-oop 2 +x x ),0(x)}. (207) 



Therefore the explicit form of the classical equations ( P05| ) for the wave created 
by the collision is 

D^Fni = 0, (208) 

.2 2 - 2 

D^F^ = <5(-x.)[-a;*pi,-oopi] a;i ViV l (a;_L)[-oopi, -x*pi] XJ _, 
s s ~ s 

2 2 - 2 

D^F.f, = 8(-x*){-x.p2 7 -oop 2 } x± \7iU t (x±){-(X)p2,-x.p2}x ± - 

s s ~ s 

These equations define the classical field created by the collision of two 
shock wavesj^j Unfortunately, it is not clear how to solve these equations.^ One 

"They are essentially equivalent to the classical equations describing the collision of two 
heavy nuclei in Ref. 55. However, we do not impose the additional boundary conditions at 
*f = 0. 

'"In Ref. 56, the numerical solution was suggested. 
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can start with the trial field which is a superposition of the two shock waves 
( |170| ) , and improve it by taking into account the interaction between the shock 
waves order by ordertil The parameter of this expansion is the commutator 
g 2 [Ui, Vk]. Actually, there are two independent commutators, 

Li = L\t a , L\ = if abc UfV bj , 

L 2 = L a 2 t a , L% = ie ik f abc U bi V ck , (209) 

where Cik is the totally antisymmetric tensor in two transverse dimensions 
(ei2 = !)• In these notations = L\ and [£/i,Vfc] — (i *-> fe) = ti k L 2 . It 

can be demonstrated that each extra commutator brings a factor In (each 
commutator means higher term in the cluster expansion in Fig. |30| ), thus this 
approach is a kind of LLA. It is convenient to choose the trial field in the 
form f\ 

4 0) = 4 0) = 0, A| 0) = 6{x.)Vi + 9(x*)Ui + 0(x m )9{x.)Ai (210) 

where Ai(x±) = Ui(x±) + Vi(x±) + Aj(xj_) is a pure gauge field satisfying the 
gauge condition 9,Aj — i[A i; Aj] = 0. The explicit form of A, is 

A'Ocl) = ig ^(U^U + V^V-^)L 2 + 0(L>) (211) 

= -ig [ dz / %~ Z j!* {U x Ut + V x Vt - l)L 2 (z x )dz x + 0(L 2 ). 
J 2n(x -z) ± 

In the first nontrivial order one gets: 

r(i) if, 1 



^ = -2^J d ^ -xf l +(x~z)l +i e A ^ ) 



g f, e ik {x- z) k ( {x-z)\ ! 



A ' ] = #2 f dz ± —!—ln(-x 2 + (x--z)i+ie)L 1 (z 1 _), 
lbn z J x* + ie 

- [ dz± 

I67H J x. + ie 



A * 1] = -t!z3 ^TT ln(-a;f + (f-z)i + ie)Li(z ± ), (212) 



x In the paper of Ref. 3, I used a slightly different trial configuration A*°' = A^ ' = 0, A^ ' = 
6(x.)Vi + 9(x*)Ui. The difference Aj is corrected by the AW term, so the results for the 
total field A(°) + AW are the same. 
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where 



„2 — 4 



longitudinal part of x 2 . These fields are obtained in the 



background-Feynman gauge. The corresponding expressions for field strength 
have the form 



F, 



(i) 



/ k 



F. 



(i) 



F 



(i) 



dz± 



e ik / dz 



— x?, + (x — z)\ + ie 
1 



gs_ i , i 

4tt 2 

9 
2tt 2 

87H / -i; + (a; — z) , + ie V x* — it 



-A2 — ^^2(^_l), 
atf + (a; — z) + 16 



i[^U (0) ]. 



(213) 



g ik Li(z±_) e ik L 2 (z±) 



a;* + ie 



(x - z) k 



-xl + (x — z)\ + it 



x, — ic 



x, + ie 



In terms of usual Feynman diagrams (when we expand in powers of source 
just like in Sec. 4.2) these expressions come from the diagrams shown in Fi g. p2| . 
When we sum up the three contributions from the diagrams in Figs. |32ja, |32|b, 





(a) 



(b) 



(c) 



(d) 



Figure 32: Pcrturbativc Feynman diagrams for the field strength (213) 



and |32| c the three-gluon vertex in Fig. |32| a is replaced by the effective Lipatov's 
vertex j30| ) and we get fl213| ) up to the terms jpdid k U k and jpdjd k V k standing 
in place of Ui and Vj. However, as we have discussed in Sec. 3, the difference 
Ui — -gzdidiJJ k = g!jp[Ui,Uk] (which has an additional power of g) will be 
canceled by the next-order perturbative diagrams of the Fig. |32] d type. 
Let us now find the effective action 



S cS = S QCD {A) + / <rx ± V m (x x )Y?(x ± ) + / d d x ± W m U?(x±) (214) 
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in the semiclassical approximation. In the trivial order the only non-zero 
field strength components are F.y = 5(^x*)Ui(x±_) and F^' = 5(-x,)Vi(x±), 
hence we get the familiar expression = J d 2 x±V al U^ . In the next order 
one has 

S (i) = J^^P^^P^^^^ 

+ 2 Jd 2 x± J du(TrV % ([-oopi, upi] x F,i(up\ + x±)[upi, -oopi] x ) (1) 

+ Trll 1 ([-oop 2 , up 2 ]xF*i(up 2 + x^)[up2,<xp2\x)^ ) ■ (215) 
Above, we have seen that the effective action contains In (see Eq. (|l94|)). 



With logarithmic accuracy, the right-hand side of Eq. (215) reduces to 



S« = -- / d i xFP a \x)Fil )a {x). (216) 



J d 2 x±_2TrLi{x^) ([x±, —oop 2 + Xj.]^ - [x±, —copi + x±] 



(i: 



The first term contains the integral over d 4 x = ^dx,dx*d 2 x±. In order to 
separate the longitudinal divergencies from the infrared divergencies in the 
transverse space we will work in the d = 2 + 2e transverse dimensions. It is 
convenient first to perform the integral over a;* determined by a residue in 
the point x* = 0. The integration over remaining light-cone variable x» then 
factorizes in the form / °°(ia:,/x, or /_ dx,/x,. This integral reflects our 
usual longitudinal logarithmic divergencies, which arise from the replacement 
of vectors n and n' in ( |193| ) by the light-like vectors p\ and p 2 . In the mo- 
mentum space this logarithmical divergency has the form J da/ a. It is clear 
that when a is close to a (or a') we can no longer approximate n by p\ (or n' 
by p 2 ). Therefore, in the leading log approximation this divergency should be 
replaced by In ^ , 

i />oo y r a ^ Q 

dx,— = I da > / da- = In—. (217) 

o x, J a J a a a' 

The (first-order) gauge links in the second term in the right-hand side of 



Eq. (216) have the logarithmic divergence of the same origin, 
r id) * f° dx * f i2 r ( e ) 

[lJ " - 0O! " +lJ = -a^J^—J dz± W^lff LlM 

^.-» + .J» - £ fj± J A^*), (218) 
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which should also be replaced by In -p- .^| Performing the remaining integration 
over x± in the first term in right-hand side of Eq. (216) we obtain the the 
first-order classical action in the form 



5 (D 



8tt 2 n cr' 



(220) 



or 



fxxtfv^LUxjLUvi) + L a 2 (x ± )L a 2 (y ± )) 
5 (i) = ^! in ° [ d 2 [ + £« j_LS ] 



r(e) 



(221) 



Note that in the trivial order the three terms in Eq. (214) are equal up to the 
different sign of the S(A) term. It can be demonstrated that this is true in the 
first order, too: 



d 2 Xj_2TrV% 



(o+i) _ 



= / 2Tr wf +1) [/ 2 = -S(A) 



A more accurate version of Eq. (221) has the form (see Appendix 7.5) 

,2 



1) 



2tt a' 
1 



d 2 xj_ 



(222) 



(223) 



1 sab, 



L\—L\ + L a 2 (U^U + V^V - —) a °L: 



dl 



L a(^_ u ^ u _ u ^ v ) L bjk 
lK dl dl J 2 



&l 



&l 



u 



V) ah L\ 



d 2 ± 



0({U,Vf). 



v The fields A. and A* in Eq. (212) look like they satisfy the condition x„A, + x.A t = 
implying the fact that P exp ig J due^ A^(un + x±) = for any vector e = cpi + CP2- One 
may suspect that the proper limit at e 2 — > is to set [x± , — oopi + Xj_] and [x±, — oop2 + *x] 
to 0. However, careful analysis with the slope of the Y operators n = api + ap2 instead of 
pi and the slope of W operators n' = cr'pi + <x'p2 instead of p2 shows that 



[x_L, — coe + a; j 



16?r 2 



d 2 zj_ 



(x - z)' 



'ct'/ct' + c/c, fa' 
■ In 



WW 



-Li(z ± ) 



f/f+cr/cr co- 
■ In 



(219) 



/ct'-c/c ?cr' c/c-ct/ct 



leading to (21£) if c — * tr or c — > c'. 
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It is easy to see that in the case of one weak and one strong source this 
expressions coincides with ( |200| ) (up to the terms of higher order in weak source 
which we neglect anyway). 

At d — 2 we have an infrared pole in which must be canceled by the 
corresponding divergency in the trajectory of the reggeized gluon. The gluon 
reggeization is not a classical effect in our approach, rather it is a quantum cor- 
rection coming from the loop corresponding to the determinant of the operator 
of second derivative of the action 

" * f S QCD + [d 2 x x V ai (x ± )Y?(x ± ) + fd 2 x x W"T';u>_ 

A=A 

(224) 

The lowest-order diagrams are shown in Fig. and the explicit form of the 



*•-. •••o w 




(a) (b) 

Figure 33: Lowest-order diagrams for gluon reggeization. 

second derivative of the Wilson-line operator is 

/OO nil 
du I dv[SAi(upi + ^7{5Ai(vpi + x±)], 
-co « — OO 

5W l = i I du dv[Ai(up 2 + x±),Vi8Ai(up2 + x±)}. (225) 



OO J — OO 



Now one easily gets the contribution of the Fig. p3^ diagrams in the form 



S T = -— r-ln— / d x^d yj_ 
57r J a' J 

x (vf(x ± )U ai (y ± ) - Vtix^ix^)) 



(226) 



T 2 (l + e) 



((f-^)i)(l+2e)' 
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A more accurate form of this equation reads: 

<? 3 2 N C a f 2 2 r 2 (l + e) 

x { - V?(x x )U ai (x x ) + ±-(v i {xx){U(xx)Ui(vx) 
+ V(xx)VHy x ) - lWiyx))^} + 0([U, V\), 

where O aa =TrO in the gluonic repr esent ation. In the case of one strong 
and one weak source it coincides with ( |l97| ) (up to the higher powers of weak 
source). 

The complete first-order (= one-log) expression for the effective action is 
the sum of S(°\ S*«, and S t , 

if CT ,2 f T(e) 



S cfi = / d 2 xV al (x)Ut(x) + /a In- \d*xd*y\ - _ v ' (228) 
J 8?H a' J ^ {x-z)f 

x (Ll{x)L\{y) + L a 2 {x)L b 2 {y)(utU y + VjV y - l) oiA 
2 e ij {x- z)i(z-y)j 



d l z 

n(x - z)\(z '- y)\ 
x (Lt{x)(UtU y -U~ V) ab L b 2 {y) - L%{x)(UtU z - U <- V) ah L\{y)) j 

+ ^(v i (xx){U(xx)UHyx)+V(xx)V^yx)-l}U i (yx)) aa \. 

I n th e case of one weak and one strong source this expression coincides with 
( [20 1| ) up to the higher powers of weak source. (As we discussed in Sec. 4.3, the 
new nontrivial terms in the case of two strong sources start from [Y,V] 3 \n 2 

As usual, in the case of scattering of white objects the logarithmic infrared 
divergence ~ - cancels. For example, for the case of one-pomeron exchange 
the relevant term in the expansion of e eff has the form 
2 /• 

— 9 — In - t \ d 2 xxd 2 yxf dam (VfU m ig ik + VfU? - V^U^x) 



16?r 2 & 



(x - y) 2 l 
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g 2 N c 



16tt 3 



In- / d 2 x ± V«(x ± )U a \x ± ) / d 2 y ± d 2 y ± (V b (y ± )-V b (y ± )) 



r 2 (i 



((y-y')i) (1+2£) 



(U b Hy ± )~U b J(y' ± )). 



(229) 



It is easy to see that the terms ~ - cancel if we project Eq. (229) onto colorless 
state in t-channel (that is, replace V al V b by V cl Vj). It is worth noting 
that in the two-gluon approximation the right-hand side of the Eq. ( 229 ) gives 
the BFKL kernel @. 

As an illustration, let us present the next-to-leading contribution to the 
effective action ~ [U, V] 3 In -p- coming from the diagrams of Fig. |34| type. 




Figure 34: Typical next-to-leading order contribution to S e g. 



ScS = 3 3 /a6 C ln^7 J dx±dy±d±z[K 1 (x±,y±, z ± )L^(x±)L b 1 (y±)L2{z±) 
+ K 2 (x ± ,y x ,zx)LUx±)L b 2 (y x )L c 2 (zx)], (230) 

where 

K i (x,y,z) = J ^^if i (p 1 ,p a ,-p 1 _p 2 ) e * 1 -(>-*)+*«-(»-«) ) 

V ( \ i £ikP\P2 fi 2 Pi , 2 Pi , 2\ 

i tikPiPl ( 1 , Pi , 1 i Pa 



K 2 (P1,P2,P3) = —7-* 2"^" ~~2 2 ln_ 2 + ~~2 2 ln ~~2 ' ( 231 ) 

Effective action as integral over Wilson lines 

In this section we will rewrite the functional integral for the effective action 
(193) in terms of Wilson-line variables. To this end, let us use the factorization 
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formula ( |189| ) n times as shown in Fig. |35|. The effective action factorizes 
then into a product of n independent functional integrals over the gluon fields 
labeled by index k: 



e iS eB (.U,v-,v) = 



JdA 1 . . . DA n+ i expz{^r„Vi + S(An +1 ) (232) 
n+l^n + S(A n ) + ... + W 2l Y{ + SiAx) + WtUi}, 



+ W, 

where the integrals over x± and summation over the color indices are implied. 
As usual, Y£ = ^d^k and W l k = ?-Wld l W k where 



Yk(xj_) = P expjg / du n^Ak^{un k + 

J — oc 
/>oo 

W k (x±) = P expig du n^Ak^un^ 1 +x ± ), (233) 

J OO 
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and the vectors n k are ordered in rapidity: rjo > r\ n > f]n-i ■ ■ ■ V2 > Vi > Va- 
To disentangle integrations over different A k we use the formula 

e iJdx ± WiY* = te b ( i a i -igW i )(ff i -igY i ) (234) 
DV(x ± )DU{x 1 _)e i I dxxWiUi+i I dx±ViYi - i I dxxViUi . 



The determinant gives the perturbative non-logarithmic corrections of the same 
order as the corrections to the factorization formula (190). In the LLA they 
can be ignored, consequently, we obtain 

e iS Ktt (u,v) = J DAl . . . DA n+1 DU 1 DV 1 . . . DU n DV n 

x exp i{ViY^ +l + S(A n+1 ) + W l n+1 U n>i - V nti Y* + ... 
+ WtiUi - V 2 l U 2l + V 2 ,iY* + S(A 2 ) + W 2 ,iUi - V hi Ul 
+ V ht Y 1 l + S(A 1 ) + WlU i y (235) 

Now we can integrate over the gluon fields Ak, 

' DA k e Vk -* Y * +s( - Ak)+Wk '* u i<-i = e ^ ri f(V fe ,c/ fc - 1 ;AT;) . ( 236 ) 

at sufficiently small A77 

S eS (V k , V k ~ x ; At?) - VkM-x - ^vK{V k , U k -±) + 0(Ar/ 2 ), (237) 
where K is the kernel calculated in the previous section, 

K(V, U) = -a s I d 2 x ± 

1 1 1 , 

2 



I%e ik (rf$-U%--U~V) ab L i 1 



1 

-in 



d\ dl 

(Vi (U^(hxd\)U + y f (In dl)V- (In U l ^j °° | . (238) 
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Performing the integrations over A k we get 



JS ott (U,V) = 



J DVxDUt . . . DV n DU n exp [iV t W n + K(V, U n )A V (239) 

- iV n ,iU l n + Wn,iUi-l + K ( V n, Un-l)A V + ...- iV 2i U l 2 

- iViUu + K{V 2 , U 1 )Ar ] - iVfUu + iVfUi + K(Vi, Z7)Ar?|. 

In the limit n — > oo we obtain the following functional integral for the effective 
action 



iS et t(U,V) 



DV(rj)DU(r]) 
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exp ■ 



{iV?U ai (r]) (240) 

u(n' )=u 

d V [ - iV ai (r,)U?( V ) + K(V( V ), U(r,j) }. 



where we displayed the color indices explicitly. This looks like the functional 
integral over the canonical coordinates U and canonical momenta V with the 
(non-local) Hamiltonian K(V, U). The rapidity ij serves as the time variable for 
this system. Let us demonstrate that perturbative expansion for the functional 
integral ( |240|) determines the effective field theory f or re ggeized gluons. To get 
the perturbative series for the functional integral (240), we write down U{rf) 
and V(rj) as 



U(x±,rj) = e - ia ^ x ^\ 



(241) 



((p a (x±,rj) and ir a (x±,ri) are scalar fields) and expand in powers of g. In the 
leading order in g we obtain 



j Dir( V )D4>(ri) 


exp < 




4>W )=<k V 



exp | -id^d^im) 



a s [dir(r]), 8<f>(ij)]-=r-[dir(r)), d(t>{rj)] 



(242) 



where d = d\ + id 2 , d = d\ — id 2 . The bare propagator for these fields is (cf. 
Ref. 2) 



(<P{x±,v) <t>(y±,n')) = °) 
(4>(xx,r)) ^{vwn')) =9( , n-v')((x±- 
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(n(x±,T)) n(y±,r/)) = 0, 
i 



PI 



y± 



(243) 



The 9 function in this formula satisfies the condition 9(0) = as can be 



easily seen from the limiting formula ( 239 ) . It is convenient to include the 
g 2 7r<9^ ln<9^0 in the kinetic term rather than in the interaction Hamiltonian. 
Since this expression is IR divergent one should at first consider the regularized 
S e fi 



a s [d<l>M^—Ad<t>M} - td^d^M \ (244) 

and then take the limit /i 2 — * 0. (Alternatively, one can use the regularization 
d = 2 + e for the number of transverse dimensions as it was done in Sec. 5.1.) 
The propagator takes the form 

(4>{x u rj) <f>(y±,v')) = 0, Hx^rj) n(y ± ,r)')) = 0, (245) 

dp± 



e^-«)x Wa;±)??)7r(y±)77 / ))=% _^ ) ^ e -ft^(W)in^ 



which coincides with the propagator of the reggeized gluon (|73|). 
Since the only non-vanishing Green functions are 

(<f>{xuri) ■ ■ ■ <f>(x m ,ri)'K(yx,r)') . . . ir(y n , ?/)) 

with m = n, the number of reggeized gluons is conserved. It is easy to see that 
the Feynman rules for the Green function 

{(t>{xx,rj) . ..(f>(x n ,rj)Tr(yi,rj') . . . Tr(y n , rj')) 

reproduce the diagrams for the quantum mechanics of n particles with Lipa- 
tov's Hamiltonian ( |75| ) (see Fig. 10). 

In the next order in the expansion ( |24l| ) we get 



j Dir{n)D<l>{n) 


exp I 




<l>{v' )=4> V 



(246) 



<t>(v) 



[90, Bit) [90, Bit] } + i^-K m (0, tt) + £-K w (0, n) | , 
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where 



#(3)(</>,7r) = \ [[d i (l),(t>},d l Tr]-T-[d j (f>,d J TT) 



(247) 



and 



([[9,0, 0],a,-7r] + 2[0, [9,0, a i7 r])^-([^, 9,-tt] - (t - j)) 



[0^,0,.*] J^J^t", [0^,0 fc 7r]] - (i «-> fc)) L + {tt <- 0} 



0* 9j_ 



#(4)(0,tt) = [[^0, 0], [^7r,7r]]-^-[9(/),a7r] 



(248) 



[[d<l>,<t>]Air]-±-[[d i <t>,<l>],d i n] + 



The number of reggeized gluons is no longer conserved, hence we get the field 
theory of reggeized gluons with Feynman diagrams shown in Fig. pq. In higher 



i / 

- 



T=T| 



Figure 36: Feynman diagrams for the field theory of reggeized gluons. 

orders we will get more compli cate d ir m (f> n vertices. 

It is intersting to compare (235) with Lipatov's effective action for reggeized 
gluonsBH In these papers the reggeon is defined as a scalar field depending 
on both transverse and longitudinal coordinates. The integration of Lipatov's 
effective action over longitudinal coordinates of the reggeons in the LLA repro- 
duces the first two (BFKL and three-pomeron) terms in the expansion (235). 
Hopefully, the integration of the Lipato v's a ction in the NLO LLA, NNLO 
LLA etc. will reproduce the expansion (235) order by order in perturbation 
theory. 



85 



5.3 Semiclassical approach to Wilson-line functional integral for the effective 
action 



Perturbation expansion ( |241| ) is relevant when the characteristic Ui and Vi 
inside the functional integral (24C) are ~ 0(1). However, we shall see below 
that at high energies the characteristic fields in this functional integral seem 
to be large, consequently the expansion ( [24 ip m ay be useless. In this case, we 
can try to calculate the functional integral fl240| ) semiclassically. The classical 
equations for the functional integral (|240D are 



(idi+gWW 

{id t +g[U t )V l 

with the initial conditions 

U(t)) = U at r? = r]' , 



:K(U,V), 



K(U,V), 



WSU 



V(r/) = V at t] = rjo. 



(249) 



(250) 



Let us denote the solution of these equation by U(x±, rf) an d V(x ±, rf). In the 
LLA the semiclassical calculation of the Wilson-line integral ( 240 ) is equivalent 
to the semiclassical calculation of the original functional integral (193). I will 
make a conjecture that the saddle point of the original functional integral (202), 
satisfying the classical equations (208), corresponds to the classical solution 
(189) of the Wilson-line integral (24C) even beyond the LLA: 



exp {iVi? 1 M + iWibQlT + iS(A)} 



(251) 



exp J iVlUito) + f° dJ - iV*(ri)Ui(ri) + K(V( V ), U(rj? 
[ J n' J 



where A is the classical solution of the equations ( 208 ). (As in previous section, 
we do not display the integrals over the t rans verse coordinates). Being a 
quantum correction, the gluon reggeization ( |227| ) exceeds the accuracy of the 
semiclassical approximation, hence we can drop the last (reggeization) term in 
the kernel dgggl). 



Talking the variational derivative of both sides of Eq. (251) with respect 
to V, we obtain 

YM = Ui{n). (252) 

If we now take the derivative of both sides with respect to 770, we get the 
equation 

iViYi^o) = iViUi(rio) = K(V, Ufa)), (253) 
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which may be used for the calculation of K. Correspondingly, one can differ- 
entiate with respect to r]' resulting in 



iV i (Vo)Ui = K(V(r 1 ' ),U). 



(254) 



Since V in Eq. (253) and U in Eq. (|254j) arc arbitrary, we may substitute V(r]) 
and U(rj) instead : 



iVi(ri)Ui{ri) = -iVi(ri)Ui{rj) = K(V( V ),U( V )). 



(255) 



The exponential of the Wilson-line functional integral vanishes except for the 
non-integral term ViU l (r] ) = ViY % (r]o), so 

exp {mY^rj'o) + iW t ( m )U l + iS(A)} = exp {iV^M} 

= expliVii-n^Ui}. (256) 

Thus, in a semiclassical approximation (and with the assumption mentioned 
above) we obtain 



S eS = ViUifro) = Vi(rQUi = -S(A), 



(257) 



so that all the three terms in left-hand side of Eq. (251) contribute equally up 
to a different sign for S(A). We have checked it in LLA and it is crucial to 
check it in the next-to- leading order. From Eqs. (255) and (257) we see that 
the effective action in the semiclassical approximation can be written down 
also as 

S cS = Vi(ri)Ui(Ti) (258) 

for arbitrary 77. 

Instead of taking variational derivatives of the kernel K(V, U), it is possible 
to calculate U = Y directly. One obtains (cf. Eq. fl218| )) 



[x±, -oopx + x± 



[oopi + xj_, xj_ 



l(i) 



l(i) - 



^ln^ 
2tt a' 

1 



((- 



2tt a 



d 2 z ± (259) 
(i 1 (z i ) + 2[C/ l (z i ),A l (z ± )]), 

d 2 z± 



ut^.u + rti{& ± u)± zy b L\{ Zl _) 
P I P I P I >' 



PI 
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and, therefore, 

[oopi + x±, — oopi + ajjl^ = — \n— d 2 

it a' J 



z± 







1 




[((- 


— o 

pi 



-a ' 



p± 



l -^t a hx^ J (fz±^([x±\-U^^-i{d k U)^-\z ± ^ ab L b 1 {z±) 

P± P± 



(260) 



The derivative W U is half of the coefficient in front of In in this formula so 
we obtain 



^U=^[m^-(d k U)\ ±-L\-^[U^U\ e ik Ll, (261) 



gk 



ab 



i 



id 2 d % ( - r t d k 



2^ \ d\ 
and, similarly, 



J_L5-^-5-(vt^.v'] gifci*, (262) 



where ?7 ee [/(??), Fee 7(77). 

For illustration, let us present a first few terms in the semiclassical expan- 
sion of the effective action, 



S eS = j d 2 Xl y,U l (263) 
2tt a' J V 52 02 



^ In ' 



2 V2tt 



((^-W(ft-W) bC {(T/t|(^F)i) C ^ + ^^A*}. 



Once we know the solution of the Wilson-line classical equations (261)- 



(262), it is possible to restore A. Suppose we want to find A(r) x , r, x±) where 



r = x?, and rj x — In 2t , Let us insert two factorization formulas at 77^ + 
and r\ x — Srj and integrate over the fields in the regions 770 > r) x + <5ry and 
Vx ^ Stj > r) > r]' semiclassically. The final integration over the region of 
rapidities rj + drj > r) > rj x — Srj takes the form 

JdA exp {iV\r, x + Ar})Yi( Vx + At?) + iW l (rj x - Ar^U^ - A v ) + iS(A)} . 

(264) 

(Here rj x + Aij denotes the argument for the classical solution V % and the 
direction of the Wilson line for Yi). Comparing this to Eq. ( |202| ), we find that 
the field A(rj x , r, x±) is given by expressions ( |212[ ) with U U(j]), V — * V(rf). 
Unfortunately, the accuracy is again up to [U(rj),V(r])] 2 . Still, we see that the 
fields contain logarithms of r\ x coming from and U (rj) V(rf) so our assumption 



about large characteristic fields in the functional integral (193) is justified. 
Note that for the infinite Wilson line in r\ x direction we can get an (almost) 
explicit expression in terms of U — > U (77) and V — > V(rj) without the restriction 
[U(ri),V(r])] <C 1. It is easy to see that 

[x± — oon v , x±+oon v ] (idi + Ai(x± + confix ± +0012^, x± — oon^] = Ai(x±, rf), 

(265) 

where A, (x±,r]) — U(x±,r)) + V(x±,r]) + A(xj_,ri) is pure gauge field satisfying 
the equation 

(idi + [Ui + Vi,) A t = 0, (266) 



(see Eq. (211). Indeed, let us try to calculate the l.h.s. of the Eq. (E6q). At 



small Srj all the contributions coming from [x± + oon,,,a:j^ — oon n ] contain Srj 



(see Eq. (219)), hence they are small. The only non-vanishing contribution 



comes from Ai(xj_ + oon^) which coincide with Ai(x±,rj) in the background- 



Fey nman gauge (266) 



6 Conclusions and outlook 

First I would like to discuss the relation of this method to other approaches to 
the high-energy QCD discussed in the literature. 

By far, the most popular approach to high-energy pQCD is the direct 
summation of Fcynman diagrams (and related methods based on unitarity 
relations in s and t channels). Although the majority of the results in pQCD, 
including the NLO BFKL kernel, were obtained by this method, I think that 
even in pQCD, the Wilson-line language, combined with the calculation of 
the propagators in the shock-wave background, is technically more powerful. 
(Perhaps the comparison of the diagrammatic calculation of the three-pomeron 
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vertex in Ref. 45 to the computation of the gluon propagator in the shock-wave 
background in Sec. 7.3 demonstrates this most clearly). 

The dipole pictureEZl has an advantage of visual interpretation of the high- 
energy scattering, especially in the case of DIS at small xcilEa The dipole 
language is a light-cone version of the Wilson-line approach combined with 
large- N c approximation for the wave functions at small x. However, it is hard 
to think about the effective action in terms of the dipoles, since in order to 
study the energy evolution of the effective action we must take into account 
not only the creation of the new dipoles, but their multiple creation and re- 
combination, which is difficult to define in the framework of the dipole model. 

The most close in spirit to our semiclassical method is the renormalization- 
group approach to the high-energy scattering from the large nuclei advocated 
in the papers of L. McLerran and collaborators (see e.g. Refs. 4, 52, 58). 
In this approach, the small-x evolution of one strong shock wave (created 
by a source p{x±)) is studied in the light-like gauge. With such a choice of 
gauge, the second shock wave can be treated perturbatively at the very end of 
the evolution process. In our terms, this amounts to the solution of classical 
Eqs. (208) using the trial configuration Aj = Ui9(x*) (instead of starting point 
Ai — Ui9(x*) + Vi8(x m ) + Ai taken in this paper). Unfortunately, due to 
different gauges adopted in our paper and Refs. 52, 58, the treatment of the 
boundary terms in the functional integral is different, leading to the different 
sources for the shock waves and making hard to compare the intermediate 
formulas. However, since the first-order (BFKL) results coincide I think these 
effective actions are essentially the same. 

In conclusion I would like to outline possible uses of this approach. The 
ultimate goal is to obtain the explicit expression for the effective action in all 
orders in In . One possible prospect is that due to the conformal invariance 
of QCD at the tree level our future result for the effective action can be formal- 
ized in terms of conformal two-dimensional theory in external two-dimensional 
"gauge fields" Vi and Ui. So far, I was not able to use the conformal invariance 
because it is not obvious how to implement it in terms of Wilson-line opera- 
tors. We can, however, expand Wilson lines back to gluons. The conformal 
properties of (reggeized) gluon amplitudes are now well studied. In the coordi- 
nate space the BFKL kernel is invariant under Mobius group and therefore the 
eigenfunctions of BFKL kernel are simply powers of coordinates. It is not clear 
which part of the conformal symmetry survives for the full effective action, yet 
there is every reason to believe that it will simplify the structure of the answer 
even after reassembling of Wilson lines. 

The semiclassical approach developed above for the small-x processes in 
perturbative QCD can be applied for studying the heavy-ion collisions. As 
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advocated in Ref. 4, the coupling constant for the heavy-ion collisions may 
be relatively small due to high density. An estimation of the corresponding 
"parton saturation scale" Q s gives ~ 1 GeV for RHIC and ^2 — 3 GeV for 
LHCP so g(Q s ) is a valid perturbative parameter. On the other hand, the 
fields produced by colliding ions are large, so that the product gA is not small, 
showing that the Wilson-line gauge factors V and U are of order of 1. Thus, 
we have a perfect situation to try sQCD methods. 

It should however be mentioned that in this paper we considered the special 
case of the collision of the two shock waves, namely without any particles in 
the final state. It follows from the usual boundary conditions for Feynman 
amplitude (81) which we calculate: no outgoing waves at t — > oo and no 
incoming fields at t — > — oo (the latter condition is satisfied automatically by 
the A\ t ^_ ca = choice of gauge). However, people are usually interested 
in the process of particle production during the collision (see e.g. Ref. 59) 
since it gives the experimental probe of quark-gluon plasma. In this case, 
our approach must be modified for the new boundary conditions — we must 
solve the classical equations (208) with Feynman boundary conditions only at 
t — > — oo. The boundary condition at t — ► oo depends on the problem under 
investigation: in the case if we are interested in the the total cross section (cut 
diagrams) we must calculate the double functional integral corresponding to 
the integration over the "+" fields to the right and the "-" fields to the left 
of the cut (see Ref. 43). (This is actually a functional- integral formalization 
of Cutkosky rules). In this case we may use the usual (Feynman and c.c. 
Feynman) propagators for each type of the fields. The boundary condition 
requires that two types of the field — the left-side "-" fields and the right-side 
"+" ones — coincide at t — > oo. (This boundary condition is responsible for 
the S(p 2 )9(po) propagators on the cut). Finally, to find the total cross section 
of the shock-wave collision in the semiclassical approximation, we must solve 
the double set of classical equations for "+" and "-" fields with the boundary 
condition that these fields coincide at infinity (cf. Ref. 60). The study is in 
progress. 
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7 Appendix 



7.1 Wilson lines from Feynman diagrams 

Let us demonstrate that the relevant operators are Wilson lines (|J) . The typical 
contribution to the Green function of the fast-moving quark (with ctk <C o~) 
is shown in Fig. |3?] where the gluons have a <C o. Consider the loop integral 



k k-p k-p' 




Figure 37: Typical diagram for the propagator of fast-moving quark. 



over p. Since we can neglect a p as compared to ctk , the quark propagator with 
the momentum k — p reduces to 

Here we have used the fact that c/ M „ in the numerator of the gluon propagator 
connecting the lines with very different rapidities (= a's) can be replaced by 

I will prove now that if I replace the propagator ( |267D by 

—(Jr. (268) 

—p p + teat 

the value of the loop integral over p remains unchanged. Indeed, the integral 
over p is the sum of the residue in the pole corresponding to the fast-quark 
propagator ( 267 ) and/or the residues in the slow-gluon propa gator s. Let us 
consider both residues in turn and verify that the replacement ( |26§| ) does not 
affect the residues. 

First, if I take the residue in the pole 

ft=A-^i (269) 
ct k s 
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corresponding tho the quark propagator, the typical slow-gluon denominator 
takes the form 



(a p + d p )(P p + P p )s - (p + p) ± 

[ptp + d p )/3 p s - [p + p)\ + (a p + d p )/3 k 



(270) 



-(k-p) 



The hrst two terms are or order of m 2 while the second two ones are ~ —to 2 
and hence they can be neglected, which corresponds to taking the residue at 
the pole P p — in the propagator ( |268| ). (Here we have used the fact that 



P>< ~ see below )- 



Second possibility corresponds to the residue taken at 

(p + p)i 



P P = -P P + 



(dp + d p )s 



(271) 



in one of the slow-gluon propagators. The quark propagator (267) then takes 
the form 

^ 2 . (272) 



Pp + (c,„+«.,)s + Pk 



iea k 



(a p +a p )s ' a k s 

2 

Again, the hrst two terms in the denominator are ~ ^— , while the second two 

° ' a p s ' 

2 2 

ones are ~ ^— <C IIL - and can be neglected which is exactly equivalent to 

OLp S OL fc s 

replacing the Eq. ( p67| ) by Eq. ( p68| ). 

Henc e, we have proved that the propagator of the fast quark can be re- 
duced to (268) which is nothing but the cikonal gauge- factor in the momentum 
representation. 



7.2 Quark propagator in a shock-wave background. 

Let us now find the quark propagator in the shock-wave background. We start 
the path-integral representation of a quark Green function in the external field 
B n , 



((■• 



V = -i 



dr\\ x 

drAf- 1 



Ve 



x{t)—x 



x(0)=y 



Vx{t){- /b+ B n (x(r))}e 



Jo 



x P exp{ lff / dt(B^(x(t))x»(t) + -a^G%(x(t))}, (273) 
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where = §(7^7^ — 7^7^)- First, it is easy to see that since in our external 
field ( p0| ) the only nonzero components of the field tensor is G n oi only the first 
two first term of the expansion of the exponent exp{Jdt^(aG )} in powers of 
(crG) survive. Indeed, a^G% = j^p%l l G^ and therefore (crG a ) 2 - 2 7 4 ) 2 = 
since fi 2 commutes with Consequently, the phase factor for the motion 
of the particle in the external field (|90| ) has the form 

Pe i9 fa dtB " < - x< - t ^ ± » ( - t } (274) 
T dt'Pe 9 It> dtB "^ t ^ i ^ g G n (x(t'))Pe 9 So dt£ v K*))*^*) 



Let us consider the case cc* > 0, y* < as shown in Fig. Similarly to the 
case of scalar propagator, we can replace the gauge factor along the actual path 
x n(t) by the gauge factor along the straight-line path shown in Fig. [li] which 
intersects the plane = at the same point (z ,z±) at which the original 



path does. The gauge factor (275) reduces to 



U n (z±) + T^f^C/^) (275) 



where the last term was obtained using the identity 

-^-U(x±) = -—[dx*[(xp^\—x*pf ) ] x G oi (—x it pf ) +x 1 _) 
OXi s J so s 



[-x^Voop^, (276) 
so 



and the factor x* (r') in Eq. (274) comes from changing of variable of integration 
from t to Similarly, the phase factor for the term in the right-hand side 

of Eq. fl273|) which contains fi n (x(r)) = j- i> 2 B^ 1 (x(t)) in front of the gauge 



factor Eq. (|273|) can be reduced to 

- aJ-[-x*P ( i ] + x ± , -00 + x x ] = -^8{x,)[U(x±) - 1]. (277) 

OX* So 

(The factor ~ (crG) is absent since it contai ns ex tr a jl 9 and p% = 0). If we now 



i nser t the expression for the phase factors (274), (277) into the path integral 
djrj) , we obtain (cf. Eq. @) 



"2 



poo /■x(t)=x »o -2 

5{x,)[U n { Xl _)-l] / drN- 1 Vx{t)e~"^ dt ~ (278) 

JO Jx{0)=v 
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dr / dr' / dz5(z*)N~ 



x(t)=x 



x{t')=z 



Vx{t) /k(j)t 



x# \T j 



x(t')= 



x(0)=y 



Vx{t)x*(T')e 



Make a shift of time variable r 1 and using Eqs. ( |9q ) and (Q9J) to perform path 
integrals in the right-hand side of Eq. ( |27§| ), it is easy to reduce the path- 
integral expression for the quark propagator in the shock-wave field ( |9l| ) to 



(4W 



4ir 2 (x — y) 2 



S{x*)[U n - l]{x±_) 



(279) 



2it 2 (x-z) 4 



i J dzS(z* 



4ir 2 (z — y) 2 ' 

(ft - f)i* 



2-k 2 {x-z) 4 



2^{z-yf 

t-i 

2n 2 {z-y) i 



(in the region x* > 0, y « < 0). The propagator in the region < 0, > 
differs from Eq. ( |279j ) by the replacement U n <-> . In addition, the 
propagator outside the shock-wave wall (at a;*,y* < or > 0) coincides 

with bare propagator, so the final answer for the quark Green function in the 
B n background can be written down as: 



I' 



V = 



2ir 2 (x - y) 4 
+ i dz5(z*) 



^^m^iKz^n-y,) 



lU*-l](z±)6(y.)6(-x.)} 2 J { J y)4 , 



(280) 



where we have used the formula 



dzS(zt, 



2tt 2 (x- z)^ 2 27T 2 (z~y) 4 



2tt 2 (x - y) 



Mx.)-0(y.)) (281) 



to separate the bare propagator. 

Now, one easily obtains the quark propagator (105) in the original field 
Eq. (HI) by making back the gauge rotation of the answer ( |280|) with matrix 

n- 1 . 
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7.3 One-loop evolution: Wilson lines in a shock-wave background. 

The convenient way to get the kernel of the evolution equation is to calculate 
the derivative of the two- Wilson-line operator with respect to the slope of the 
supporting line. Formally one obtains: 

^Tr{U(x ± )UH yi _)} (282) 

= igC [udu(Tr{[oo,u] x F*.(up c + x±)[u, -oo] x U" t (yx)} 

— Tr{U(xj_)igC udu[—oc,u] y F#,(upt + yj^)[u,oc] v }\. 



The kernel is the result of the calculation of the right-hand side of Eq. (282) 
in the shock-wave background. 

Consider the operators £A and U'^ in the external field formed by slow 

gluons with a <C J Making the rescaling ( p^ ) we obtain: 

( [<X>PA , - OOPA\ x [~COPA , OOP A\ y ) A 

= ([°°Pa \-°°Pa ]x[-ooPa ) '°°P^ ) ]j/)s ; (283) 



where the shock- wave field is given by Eqs. ( p8|) - (|90|). Equation ( p82|) reduces 
to 

C^(U(x ± )U\y ± )) A (284) 

2 p 

= ig— udu([oop^\up { ^]xF r .(up^ + x±)[up^\-02p { ^} x Ud(y ± )) B 
so J 

2 t. 

- ig— udu(U(x ± )[-oop i 2 ) ,up { 2 ) ] y F^(up { 2 ) + V±)[up < A ,ocp { 2\) B - 
so J 

Since the (F* ) component of the field strength tensor (^) vanishes for the 
shock-wave field, the only nonzero contribution comes from the diagrams with 
quantum gluons. In the lowest nontrivial order in a s there are three diagrams 
shown in Fig. |3^. 

Consider at first the diagram shown in Fig. |38| a (which corresponds to the 
case x* > 0, < 0). The relevant contribution to the right-hand side of 
Eq. (|||) is: 

(°) „,J°)l *ar„.J0) (0)1 „ / , , r (0) (0)i ,b 



g / du[ocpy,upy] x t a [upy,-ocpy] x ® / dv[-oopy,vpy] y t 
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a 



b 



c 



Figure 38: Path integrals describing one-loop diagrams for Wilson-line operators in the 
shock-wave field background. 



r (0) (Q)i (( (0) . r, (0) 

x I v Pa ,°°Pa \y\\ u PA + x -t u P*{[Pm 



1 



(D a G aX 



P-P2 
Pip 



P2X r> B ri oir' P2P\ 

P ■ P2 P-P2 V z g pri + 2iF f 

X (p ( S,--p S -Po)}-v{...} Pt 
' P ■ P2 



I"! 

(0) 

V P A +^-L 



P2\ 



P-P2 P-P2 



D a G a p 



(285) 



As we discussed in Sec. 4, terms in parentheses proportional to V a vanish after 
integration by parts (see. Eq. (123)). Further, it is easy to check that since the 
only nonzero component of field strength tensor for the shock wave is G a ±_ the 
expression in br aces in Eq. (28J ) can be reduced to 00 where the operator M „ 
is given by Eq. (308). Starting from this point, it is convenient to perform the 
calculation in the background of the rotated field B n ( |9l| ) which is everywhere 
except the shock-wave wall. (We shall make the rotation back to field B in 
the final answer). The gauge factors [oo, u]£ a [it, — oo] and [oo, i>]i b [v, — oo] in 
Eq. (|285|) reduce to i a [oo, — oo] (&t b [~ oo, oo] (at > 0, < 0) and we obtain: 



g 2 t a U n ®t b U^ n du dv{ 



Xj_ 



P*OH 



vpf+y±)) (286) 

// ab 



where we have used the fact that the operator commutes with O n . Let 
us now derive the formula for the ( QO ) component of the gluon propagator 

Ux O n yjj in the shock-wave background. The path-integral representation 
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of 



I' 



y )) has the form 



go 
2p, 



1 fnarrn s 



_ (r\as~<i 

J>2^ °J)2 U $72 <p2^ ao 2 P * 



(287) 



dr I dt 
ii Jo 



x G" e 



Q n iT"V 2 ls ]JO!qQ. 



2p* 



n ao e-" v2 }\y) 



dT » e *(r'-T")V 2 



drN- 



(t)=x 
x(0)=y 



Vx{t)e'^o dt ^{^ [dr 1 f dr" 



x p e i 9/> tB ^^(*))*M(t) 5G n^( r /^p e i9/;',rfts^(^(t))iM(t) r Ar ii 

J 

x Pe 19 I^^^M^^gG^ 1 . (x{r"))Pe i9 S^' dtB °^ t ^ i ^ 



i / dr 



'Pe* 9 r' dtB " (g(t))i; '' (t) S ° gD a C Q (x(r'))Pe 9 So dts ?( ;E ( t )) i ^( t )} 



As we discussed above, the transition through the shock wave occurs in a short 
time ~ j so the gluon has no time to deviate in the transverse directions and 



therefore the gauge factors in Eq. (287) can be approximated by segments of 

(288) 



Wilson lines. One obtains then (cf. Eq. (|273|)) 



O 



-si I dr I dr' I dzSiz^Af- 1 



x(t)—x 



Vx(t)s' 



x -±-^{2[GG] n (z ± )-z[DG] i \z ± )}M- 1 f X(T ) * Vx{t)e~ l S>^ , 

x * \ T ) Jx(0)=y 

where [GG] and [DG] n are the notations for the gauge factors (128) calculated 
for the background field B^, 



[DG] n (x±) 
[GG] n (x ± ) 



du[oop 1 ,upi] x D a G^ (upi + a;j_)[wpi,-oopi] x , 
du I dvQ(u — v)[oopi, upi] x G^ n (upi + xj_) 



x [upx,vpx] x G^ (vpi +x±)[vpi,-oopx] 



(289) 
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As we noted in Sec. 4, the gauge factor — i[DG] + 2[GG] in braces in Eq. (287) is 
in fact the total derivative of U with respect to translations in the perpendicular 
directions so we get 



V)) = 2 S ° 



x Af- 1 
x Af- 1 



dr J dr' j dz5(z#) 

(r')=z X*[t') 

x(t')=Z - t 2 

Vx(t)e~ l J-r' dt -. 

(0)=y 



(290) 



Using now the path-integral representation for bare propagator ( B5| ) and the 
following formula 



drAf- 



x(r)—x | 

•/V.:/:— -r 

(0)=y ».(0) 



Jo 



. \n(x-y) 2 
167r 2 (x — y) t 



(291) 



we finally obtain the ( QO ) component of the gluon propagator in the shock- wave 
background in the form: 



w on 



dz5{z Sf ) 



\n{x — z) 2 
16n 2 x* 



(292) 



x {8iU n (z ± )Q(x*)Q(-y*) - dlU^(z ± )e(-x*)e(y*)} 



1 



4ir 2 (z - y) 2 



where we have added the similar term corresponding to the case a;* < 0, > 0. 
We need also the g|- derivative of this propagator (see Eq. ( |286| )) which is 



(x - y) 2 



(293) 



x {diU n (z ± )0(x«)0(-y*) - $lui% ± )Q(- x .)Q(y m )]-J—. 

(z - y) 2 



Substituting now the Eq. (293) into Eq. ( |28q ) one obtains 



((- 



4tt 

,-,2 



4^ 



—a 2 jT n — 



y± 



t a U n {x ± )®t b U^\ Vl _) 



a b 



i-^f/tni. ^nn u n {x±)t a ( ^ u tn {y±)t b_ (294) 
Pj_ Pj_ Q b 



which agrees with Eq. (132). 
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Let us consider now the diagram shown in Fig. p8p . The calculation is 
very similar to the case of Fig. |38| a diagram considered above so we shall only 
briefly outline the calculation. One starts with the corresponding contribution 
to the right-hand side of Eq. ( 284 ) which has the form (cf. fl285|) : 



g 2 C, I du I dvQ(u — v)[oop A y> + xx, up^ 1 + Xxl^lup^' + Xx, vp^A + 



(o) 



10) 



(0) 

x±,-oopy 



x t»[vp% 

x (jW^ + xx up*{(p M 

P2X 



xx] <8> W{yx) 

(0) <T) Vit 



-)[ 



P • P2 n V 2 g iv + 2iG^ V 2 g iX + 2iG% 



P ■ P2 P ■ P2 



1 



V 2 g pv + 2iG 



Ml 



P-P2 

P2 V 
P-P2 



P-P2 



'^V )}- V {...}p, vp^+xx} .(295) 



As we demonstrated in Sec. 4, the terms in parentheses proportional to V van- 
ish and after that the operator in braces reduce to O qo ■ Again, it is convenient 
to make a gauge transformation to the rotated field (^) which is every- 
where except the shock wave. Then the gauge factor [oo, u]t a [u, -co] in 
Eq. (295) simplifies to t a [oo, — oo]t b (at x* > 0, < 0) and we obtain 



g 2 t a U n t b ®U^ n du dv(u-v)((upf 



■Xx 



(0) 

vpy 



Xx 



(296) 



Using the expression (293) for the gluon propagator in the shock- wave back- 



ground we can reduce Eq. (E96J) to 



9_ 

An' 



t a U n {xx)t b ®U^{yx)((xx 4(<9104 
VV p ± p ± 



xx 



ab 



(297) 



The contribution of the diagram in Fig. p8| b differs from Eq. (297) only in 
change U «-> U* , x <-» y. Combining these expressions, one obtains the answer 
in the rotated field (|9l|) in the form 



16tt 3 



dz ± | [{ C/ tn (z x ) U n (xx )}* { U n (zx ) C/ tn (yx) } j 

+ {U n (xx)U^\zx)} l l {U^(yx)U n (zx)}^ 
- 6%{U n (xx)tf n (yx)}i - 5t{U^(yx)U n (xx)} 



(298) 



(x-z,y- z)_ 



(x-z)\(y- z)\ 



{U n (zx)})TT{U n (xx)U^(zx)} - N c {U n (xx)}W n (yx)\ 



{x - z)\ 
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{U n (x ± )})[U^(z ± )fTT{U n (z ± )U^(y ± )} - N c {U^(y ± )} k )1 ' 



Now we must perform the gauge rotation back to the "original" field B^. The 
answer is especially simple if we consider the evolution of the gauge-invariant 
operator such as Tr{U(x±)[x±, y±]-U* (y±)[y±, x±] + } where the Wilson lines 
are connected by gauge segments at the infinity. We have then 

(^{Tr{u c (x±-)[x±,yx}-U K (yx)[y±,x±} + }) A = 



as_ f, [x-y)\ 
^J aZ± (x-z)l(z-y)l 



x \Tr{U(x±)[x±,z±]-U^(z±)[z±,x±]+} 
x Tr {U(z x ) [z ± ,y x \- (y x )[yx,z x }+} 

- N c Tr{U(x±)[x ± ,y ± ]-U^(y ± )\y ± ,x ± ]+}\, (299) 

where we have replaced the end gauge factors like tt(oopi + x±)Cl' (oop% + y±) 
and fi(— oopi + x±)W (—oopi + y±) by segments of gauge line [x±, y±]+ and 

, — , respectively. Since the background field is a p ure gauge outside 
the shock wave the specific form of the contour in Eq. ( |299|) does not matter 
as long as it has the same initial and final points. Finally, note that the gauge 
factors in the right-hand side of Eq. (299) pr eserv e their form after rescaling 
back to the field so we reproduce the Eq. (137). 

In the general case, the evolution of the 2n-line operators such as 
Tr{UW}Tr{UW} ...Tr{UU>} come from either self-interaction diagrams or 
from the pair- interactions ones (see Fig. ^) . These pair- wise kernels have the 
form (U x = U(x±), etc.) 

c|™;p;>? = iIjb j^ ( f:™-_% ± poo 



'd( J 16tt 3 7 {x-z)\(y-z)l 

{U z }]{U y UtU x } k 3 + {U x UtU y }]{U z } k 3 - {UJ\{U V }) - {Uy}\{U x } 
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Figure 39: Typical diagrams for the one-loop evolution of the n-line operator. 



_L 



(f-f)2_( ? 7-f)2_ 



for the pair-interaction diagrams in Fig. p9| a and 

C^{E&j = -^- 3 Jdz ± [UiTr{U z Ul}-N c Ui]j^T, (301) 
for the self-interaction diagrams of Fig. |39| b type. 
7.4 Gluon propagator in the axial gauge. 

Our aim here is to derive the expression for the gluon propagator in the external 
field in the axial gauge. The propagator of the "quantum" gauge field A q in the 
external "classical" field A cl in the axial gauge e^A^ — can be represented 
as the following functional integral: 



10^0 



G%(x,y) = limiV- 1 DAAf{x)Af{y) (302) 



i J dzTr{Al(z)(D 2 g^-D"Dl 3 ~2iFf-^e"el 3 )Al(z)} 



x e 
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where = d^ — igA^. Hereafter we shall omit the label "cl" from the external 
field. This propagator can be formally written down as 



iGf v {x,y)=^ 



1 



ah 



i y)) » (303) 

where □ A " y = V 2 g p " + 2iF^ . It is easy to check that the operator in right-hand 
side of Eq. (303) satisfies the recursion formula 



1 



(Ve) 



1 



1 e 71 



.(D x F Xa —-V a ^D x F x ^) 



Ve 



(304) 



which gives the propagator as an expansion in powers of the operator D\F^ a = 
—g^t a ^ a 4>. We shall see below that in the leading logarithmic approximation 
we need the terms not higher than the first nontrivial order in this operator. 
With this accuracy 



1 



1 



-e^e" 



Ve 



{Ve) 



-V^D a F a ^ 



D x F Xri - 
ef \ 1 



We take now w 
in the form 



er e'' 

V~e + Ve 

e° 

Ve Ve)UP° x ~ v Ve 

0, obtaining the propagator in external field in axial gauge 



-D x F Xp 
-V v ). (305) 



s« 1 



1 



D x F Xv - 



Ve 



Ve 



Ve 

DxF Xp - —V p D a F 
Ve 



af3 



Ve* 



1 e a 



(306) 



where the dots stand for the terms of second (and higher) order in D x F\ p . It 
can be demonstrated that for our purposes a first few terms of the expansion 
of operators ^ in powers of Ft„ are enough, namely 



l Gf u {x,y) = {5l-V, — ) 



V 2 



1 

V 2 ' 



-2i^F^ 2 



1 



o 6v 



(S^-—V u )+... (307) 
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where the operator O stands for 



J)2 ± V-p2 ± S l ''p2 



(308) 



r P-P2 P-P2 



7" D 



7.5 First-order effective action. 

As we discussed in Sec. 5, in order to calculate the effective action semi- 
classically we can start with the trial c onfig uration fl210p . Making the shift 
A — > A + A^ in the functional integral ( |202|) , we obtain 



DA expJ / dx±Vf l (x ± )U at {x ± ) + 2 / j_ A" (a; j_ ) A ai (0, ;e_l) 



2Tr / dx± 
1 



1, 



2 [V\A t }Y 1 + (-L 1 + -[V\A l ])Y 2 



i y (D V - 2igF llv + g 2 G^) ab A bfl + 0(A 3 



(309) 



where 

Yi(x±) = [x ± + oopi,xj.] (1) , Y 2 (x±) = [x ± ,x± - oopi] (1) , 

W 1 (x ± ) = [x ± +^p 2 ,x ± }^, Wn(x±) = [x x ,x ± -oop2]^, (310) 

and the operator Q^ v is the second variational derivative of the source term 
with respect to A^, A v . The non-zero components of are 

G.. =S(-x*)(d l -i[V l ,)U^, g** = S(-x.)(d i -i[U i ,)V i ^, (311) 

while all other components vanish. In the first order in our cluster expansion 
we obtain 



c(i) _ _ 2 



((0,A? 



1 



+ 



2 5 2 



D 2 g ik - 2igF lk 



((0,L 1 ^ + ((0,[U i ,A i ]^--(( 



|0,A fe ) 



0,ii 



(312) 



El 

/3 
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\D 2 g t _ lu - ligF^u + g 2 G fl „ 



ab 



1>2 



0,Li 



P2 

a 



0, [U\A % 



Pi 



j3-ie 



0,ii 



Pi 





where i = (similarly fori) and ,A 2 | = f dz± 0,z ± jA 2 (z ± ) 

etc. We will now demonstrate that with 0[U, V] 2 accuracy one can reduce 



D 2 g l , v -2igF^+g 2 g, ll 



in right-hand side of Eq. (pi 



to ^ . Indeed, 



1 



D 2 g^ - 2igF f _ ll/ 



(313) 



f)2 + M 9 ^) 2 V f)2 



2 1 p 1 p 1 2 1 1 



It is easy to note that the term ~ does not contribute to right 



hand side of Eq. (312) because the relevant components of F^ v vanish: Fjk = 
F*. = 0. Let us prove that the last term in the right-hand side of Eq. (313) 
leads to the cont ributions ~ [U, V] 3 . Consider the first term in the right-hand 
side of Eq. (311). The corresponding contribution is -^Fit-^F^k^ + (•<->• 
*). Because F*i = Ui + O(Aj), Fj. = Vi + O(A^) this term is actually 
proportional to Ai ?7fc Afe ~ [C/,y] 3 . Let us now turn our attention to 
the second term in the right-hand side of Eq. ( [312] ). The relevant contributions 
have the structure L\ (^F.^F.ij^ - jpG^) Li^, Li^F.i^F^^Li, 
\Vi,Ai] (^F^F.^ + j^G,.^) WuAi], [V i ,A i \^F. i ^F^[U i ,A i ], 
and similar expressions with U <-» V, *<->•. All of them are clearly ~ [U, V} 3 
except the first term which is 



9 2 0,Li 



1 



13 + ie 



±F±F±~±C ± 
jyi * l ]j2 ' l jyi £)2 jyx 



1 



0,Li 



(314) 



If we neglect the [U, V] 3 terms i n clu ster expansion, the Green function in 
braces in right-hand side of Eq. (314) should be taken in the Ui9(x*) back- 
ground. This Green function has the form 



^ £)2 Y)2 F * 1 £)2 



D 2 S " D 2 



a. 



(315) 



plus the similar term ~ 6(—x*)6(y*). It is easy to see that the terms ~ 
0(x*)O(—y*) or ~ 6(— x*)6(y*) do not contribute to Eq. (314) — recall that 
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this term comes from the contraction of LiW^x) and ^1^2(2/) where both 
x*,y, < 0. 

Thus, the [U, V] 2 term in cluster expansion of Eq. (|3 1 2|) reduces to 



? (i) 
'cir 



0,A, 



0,A, 



0,Lx 



D 2 
1 1 



1 



1 1 



1 



a + it D 2 (3 — ie (3 + it D 2 a — it 
111 " 2 



s 



/3 + D 2 a 
1 1 



0, [t^A,])) +^((0,^ 



0,Li 
1 11 



1 



a£> 2 /3-ie 



0,Li 



■(( 



(Ul-.A' 



a + ie D 2 /3 
1 1 



0, fV" 
1 



(3D 2 a -it 



(316) 



0,ii 



It is easy to see that the remaining Green function connect points belonging to 
the different boundaries of the same sector in Fig. 40. It may be demonstrated 




Figure 40: Trial field configuration. 

that up to [U, V] accuracy the only effect of the background field on the Green 
function with the arguments belonging to the same sector is the corresponding 



gauge factor: (^x yjj = Qt(xx)((a 
We obtain 



w(yx), where is U, V, or A. 



([0,x ± 



1 1 



1 



a + it D 2 (3 — it 
1 1 1 



o,y± 



0,x± 



1 



-1 



1 



a + it D 2 (3 + it 



a + it p 2 + it {3 — it 
0,y x )) =Ul(p,xx ' L ' 



a + it p 2 + it (3 + it 



0,y±JJ, (317) 

0,V±ljUy, 
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((0,a? x 



1 1 1 



0,yx)) =Fj((0, 



a — ie Z? 2 /3 — ie 
In the leading log approximation Q 



1 



-1 1 



a — ie p 2 + ie (3 — ie 



o,y±}v y . 



(o, 



X± 



and 

so we get 



1 



p 2 + ie 
2/s 1 



1 (7 

— In— 6 (xx - y±) 

47T (7 



(319) 



a ± ie p 2 + ie f3 — ie 



0,y±)) = [[0,x± 



-2/s 1 



i 

2^ 



a ± ie p 2 + ie f3 + ie 

h v(H^h))' 



o,y± y 

(320) 



S. 



cir 



2tt ct' 



x <^ L 



dx ± A al (x ± )A1(x ± ) + g 2 I dx ± dy ± 



\\ Z p^ // 

^(^)((^|^(^ t ^^ + ^)|y4 ah [ c/ - A4 ] 6 ^^ 



(321) 



Finally, the effective actuion in the [U, V] 2 order in the cluster expansion has 
the form 



S, 



'$ = -^ln^2Tr| J dx±^A i {x±)A i (x±) + J dx^dy^L^xx) 
{ x ± | "4- 1 Vx]j Lx (y±) + 2L X (x ± ) ((V 1 4- 1 y±) [Pi - V 4 , A*] (y±) } 1 (322) 



((o,^|^|o )2 ,x)) 



z This formula may obviously seem confusing since 

which does not have any In ^7. However, careful analysis with the slope of the Y operators 
n = crpi + ap2 instead of pi and the slope of W operators n' = a' pi + a'p2 instead of p2, 
yields logarithmic contribution of the form 

af3 1 



(a + - ie){fl + 2-ra + ie) p 2 + ie 



0,yj 



-J-l n f-S 2 (x ± -y ± ). (318) 
4-7T a' 
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which coincides with Eq. 
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